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ABSTRACT 



The Dual Frequency Pump Method of acoustically determining 
point by point bubble cloud densities was studied to determine the 
practicality of using this acoustic technique to determine bubble den- 
sities in surface ship wakes. The dual-frequency technique of acousti- 
cally detecting bubbles utilizes a high- and low-frequency sound field 
to insonify the target bubbles. The bubbles themselves then radiate 
sound at the sidebands of the higher frequency. The frequency of the 
return sound is proportional to the bubble sizes present. The Dual 
Frequency Pump Method of bubble detection can differentiate and 
count many different-sized bubbles and is, therefore, well suited for 
determining ship wake bubble density distributions. The theory, con- 
siderations, experimental results, and recommendations of this thesis 
support the application of the dual-frequency acoustic technique to 
the ship wake problem. 



THESIS DISCLAIMER 



The computer programs provided in the appendices of this thesis 
are solely to provide insight into the possible applications of the 
equations presented within these pages. The programs are not 
verified or documented and use of these programs is at the user’s own 
risk. 
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I. INTRODUCTION 



New torpedo technology utilizes the surface ship wake and the 
near-surface region of the ocean to give their torpedos advanced 
capabilities. To better understand and counter these capabilities, the 
Naval Sea Systems Command Surface Ship Torpedo Defense (SSTD) 
program has identified the need to study surface ship wake 
characteristics (CNO project 0779) in support of U.S. weapon system 
development (Ref. 1]. In addition, both the Pacific and Atlantic Fleet 
Anti-Submarine Warfare Improvement Programs have listed surface 
ship wake data collection as a high priority action item [Ref. 2]. 

The purpose of this thesis is to investigate a new and more accu- 
rate acoustic technique of detecting small bubbles in bubble clouds 
which could be used to determine the bubble density distributions in 
surface ship wakes. This new method of acoustic bubble detection 
uses a dual frequency resonance technique originally investigated by P. 
M. Shankar and V. L. Newhouse of Drexel University [Ref. 3]. Acoustic 
bubble cloud density measurements are useful as they could serve as a 
practical verification for new optical measurement techniques being 
developed by the Naval Coastal Systems Center (NCSC) Panama City, 
Florida, in support of the SSTD program. Accurate bubble cloud 
density measurements would also aid in the development of ship wake 
computer models and provide insights into acoustic weapon 
performance in the wake region. 
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A. BACKGROUND 

The need to determine surface ship wake characteristics is not a 
new problem in the areas of Anti-Submarine Warfare (ASW) and Sur- 
face Ship Torpedo Defense (SSTD). Ever since the development of 
the torpedo, studies have been conducted on how wake characteristics 
could be used to enhance defensive or offensive opportunities for both 
targets and weapons. During World War 11, acoustic means of detect- 
ing bubbles in surface ship wakes were developed specifically for sub- 
marines and their torpedos. These techniques were crude and actual 
instances of attempting bubble density measurements were rare [Ref. 
4]. However, advances during the past ten years in acoustic-type tor- 
pedo technology have drastically increased the need to understand 
surface ship wake mechanisms and bubble density distributions. 

In this thesis, the specific area of interest with regard to the sur- 
face ship wake study is the bubble density distributions in the near- 
surface region, both inside and outside the wake. Acoustically, the 
bubble density distribution is extremely important because of the 
direct relation to sound propagation and absorption in sea water. The 
number of bubbles within a specified volume of water, that is, the 
bubble density, will determine how much sound is scattered and how 
much sound passes through that volume. If bubbles present within the 
specified volume of water have the same resonance frequency as that 
of the passing sound wave, the bubbles resonate and the effects of 
attenuation are very pronounced. If point-by-point bubble density 
distributions can be determined, then the performance of the new 
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acoustic weapons can be predicted and modified for operations in the 
near-surface regions. Measuring the bubble density distribution as a 
function of bubble size (bubble size and resonance frequency are 
directly related) and knowing the bubble density for different points in 
the wake are also critical to validating computer wake models. 

Several methods of acoustic bubble density distribution measure- 
ment have been tried in the past or are currently underway. The Naval 
Research Laboratoiy (NRL) in Washington, D.C., the Naval Postgraduate 
School (NPS) in Monterey, the David Taylor Research and Develop- 
ment center (DTRDC) in Carter Rock, and the Naval Coastal Systems 
Center (NCSC) in Panama City are just a few activities which have 
experimented with, or considered, the bubble density measurement 
problem for the Navy. 

The Naval Research Laboratory conducted a program of ship wake 
research in the late 1960s. The most noteworthy study was done with 
reverberation chambers which were used to probe ship wakes and 
determine sound absorption profiles. This research included ship 
wake geometry, wake dissolution, bubble rise, and sound absorption 
measurements. [Ref. 5] 

The microbubble measurement effort at the Naval Postgraduate 
School has been led by Dr. Herman Medwin. He is well published on 
the methods of Resonant Scattering and Second Harmonic Generation 
for bubble detection. Both methods are practical when counting bub- 
bles of approximately the same size but are somewhat inaccurate when 
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measuring bubble clouds which contain many bubbles of different sizes. 
[Refs. 6, 7, 8] 

•The David Taylor Research and Development Center was recently 
contacted by NCSC Panama City to develop an acoustic bubble detector 
to measure actual distributions in ship wakes and ship wake models. 
DTRDC proposed developing the detector in fiscal year 1988, but due 
to the lack of funds current plans for this project are now on hold. 
[Ref. 9] 

Finally, the Naval Coastal ^sterns Center has the most extensive 
ship wake research effort to date. The NCSC effort is based around a 
wake map vehicle which can make actual sound pressure level mea- 
surements looking upward through the wake to the surface. This map 
vehicle provides a profile of the wake and data that can be used to 
calculate the vertical transmission loss from the vehicle to the surface. 
However, the map vehicle does not provide horizontal wake data or 
the point-by-point bubble density distribution measurements needed 
for computer models of sound propagation. For this reason, NCSC is 
developing an optical measurement device which should be able to 
provide the point-by-point bubble density wake measurements. Still, 
this exact method of measurement has never been tried. Therefore, 
an accurate acoustic measurement device which will count all bubble 
sizes is highly desirable for verification of the optical technique. [Refs. 
10 , 11 ] 

Acoustic bubble density measurements have been made in 
association with cavitation, diver decompression sickness, contrast 
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echocardiology, and pressure flow gradients [Ref. 3]. These methods 
have been successful when measuring small numbers of bubbles of a 
single particular size. However, success has been limited when the 
bubble cloud contains many bubbles of different sizes such as those 
bubbles found in surface ship wakes. Previously attempted methods, 
which included Resonant Scattering, Doppler, and Second Harmonic 
Generation, are at a disadvantage when used for measuring a bubble 
cloud. In Resonant Scattering, a particular bubble provides peak echos 
at its resonance frequency. However, a larger bubble may reflect an 
even greater amount of energy at that particular frequency due to its 
large cross-section, resulting in a false bubble count. The Doppler 
technique by itself provides no way of separating the returns from 
different size bubbles. Second Harmonic Generation is inaccurate in a 
bubble cloud because the second harmonic radiation from large 
bubbles suffers resonant absorption by the smaller bubbles. None of 
these methods seems to provide the accurate acoustic bubble cloud 
measurements required for the surface ship wake problem. [Ref. 12] 

B. SCOPE 

A method to acoustically measure bubble cloud densities using a 
“Double Frequency Pump” sum and difference technique shows con- 
siderable potential in obtaining a high degree of accuracy when 
counting bubbles of various radii. This method involves insonifing the 
bubble cloud with two different sound fields— an “imaging” field of fre- 
quency coi and a “pump” field of frequency cop. The imaging field 
employs high-frequency sound to image the bubbles in the cloud. 
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Because of its short wavelength, the high frequency (coi) sound field 
provides good spatial resolution. The pump sound field sweeps 
through a range of lower frequencies (cop) to excite resonant 
oscillations in the bubbles. The pump frequency range is selected to 
encompass the expected range of bubble resonance frequencies. A 
bubble which is excited near resonance by the pump field in the 
presence of the imaging field will undergo large amplitude, nonlinear 
oscillations. These oscillations cause the bubble to radiate energy at 
several different frequencies, specifically at cop and coi ± cop. The 
radiated energy can be displayed on a spectrum anal 3 ^er and the 
bubble size distribution inferred from the sideband frequencies 
(cOi ± (Op) and their peak pressure amplitudes. (Ref. 121 

The Double Frequency Pump method appears to solve the pre- 
vious problems related to the inability to count bubbles of different 
sizes in a bubble cloud. It is this method that we will investigate for 
use in a simple acoustic bubble density measurement device. 

C. GOALS 

This thesis will concentrate on determining the practicality of the 
Double Frequency Pump technique for near -surface, open-ocean 
bubble density measurements. Specifically, the feasibility of using this 
technique to measure bubble densities in surface ship wakes will be 
considered. To obtain that goal, it will be necessary to complete 
several tasks. First, re-create the Double Frequency Pump experiment 
in sea water utilizing a rigorous bubble generator to generate bubbles 
with radii similar to those found in ship wakes. Part of this task will 
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require calculating the resonance frequencies of these bubbles so that 
the pump sound field can resonate all bubbles present. Second, a 
partial validation of the experimental accuracy of bubble size and 
bubble cloud density will be attempted using photographic methods. 
Third, the effects of the bubbles screening each other from the 
receiving transducer so that the energy radiated from the bubbles is 
never received need to be studied. These topics include the 
consideration of transducer types, transducer frequencies, near and 
far field beam patterns, and transducer placement for obtaining the 
best bubble count. Lastly, it will be necessary to briefly examine the 
transducer inputs required for complete isonification of the bubble 
cloud and the ideal volume in which to measure, or sample, the bubble 
density. These goals are only a few of those necessary to begin 
designing an acoustic bubble density measurement device, but they are 
enough to help show the concept practical for conducting bubble 
density measurements in the ship wake environment. 
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n. THEORY 



In order to appreciate and understand the Dual Frequency Pump 
technique for the measurement of bubble cloud densities, it is impor- 
tant to examine the theories behind bubble resonance, dual frequency 
bubble excitation, and transducer beam patterns. These are the three 
major concepts needed for this acoustic bubble detection method. 

A. BUBBLE RESONANCE 

When a sound wave strikes a bubble, the bubble undergoes radial 
oscillations. The response of the bubble depends on its size as well as 
the frequency and the pressure amplitude of the incident sound wave. 
These radial oscillations, in the form of compressions and expansions, 
are analogous to the motion of a simple, damped harmonic oscillator 
and, in fact, can be described by similar differential equations. Just as 
the damped harmonic oscillator has a resonance frequency where the 
displacement becomes maximum, the bubble also has a resonance fre- 
quency where amplitude of the radial oscillation becomes maximum. 
When a bubble oscillates at resonance, the maximum amount of energy 
is extracted from the incident sound wave. A large portion of this 
energy is then re -radiated by the bubble in all directions with the rest 
converted to heat. It is this scattering and absorption of the sound 
energy that causes large amounts of attenuation as sound passes 
through a bubbly medium. (Ref. 13] 
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As previously mentioned, a gaseous bubble in seawater behaves 
similarly to a damped harmonic oscillator, a schematic of which is 
shown in Figure 1. The differential equation for this model can be 
written by setting the sum of the forces exerted on the mass equal to 
the mass times its acceleration and then rearranging the terms. This 
well-known equation, shown below, will be used to derive an expres- 
sion for the angular resonance frequency cod of the damped oscillator. 
For the damped harmonic oscillator 

h2v H X 

m^+Rm^+sx = 0, 

where: 

m = mass 

Rm = mechanical resistance 
s = stiffness of spring 
X = displacement 
t = time. 

Dividing through by the mass m, and letting 



coo = Vs/m = undamped angular resonance frequency 



yields 



d^x Rm dx o ^ 



( 1 ) 



This equation can now be solved using the standard complex 
exponential method [Ref. 14]. The solution yields the equation for the 
natural (resonance) angular frequency of the damped oscillator. 
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Figure 1 



Simple Harmonic Oscillator with Damping 

Source: Ref. 14:7 
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cod = = damped angular resonance frequency, 

where 

X I^iii 

P = dissipation or damping coefficient. 



Several assumptions which apply to both the damped harmonic 
oscillator and the pulsating bubble now become important. First, the 
amplitude of the oscillation for both systems is considered to be small. 
Secondly, when considering a system at its resonance frequency, the 
effects of damping are small and can, therefore, be ignored. In other 
words, the dissipation coefficient is much smaller than that of the 
natural frequency. In this case. 



P = 



1 Rm 

2 m 



-> 0 



and 



cod “ cOo = Vs/m. (2) 

The damped resonance frequency is approximately the same as the 
undamped resonance frequency. [Ref. 14] 

The bubble may now be considered in terms of volume pulsations, 
V. A volume pulsation is an expansion and contraction of the bubble 
which results in a constantly changing volume. This volume pulsation, 
or so-called radial mode, is the simplest mode of bubble oscillation, 
and the one that causes the most energy radiation. Because of this, it 
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is the only mode of interest for acoustic bubble detection. Other 

modes involve shape, not volume, changes in the bubble and, therefore, 

do not radiate sufficient energy for easy detection. [Ref. 15] 

Volume pulsations of the bubble sire analogous to the displacement 

of the mass in the damped harmonic oscillator. Consider the bubble in 

Figure 2, where 

V = volume pulsation 
Vo = equilibrium volume 
V(t) = instantaneous volume 
s = bubble stiffness 
m = inertial constant 
P = dissipation (damping) 

R = instantaenous bubble radius 
Ro = average bubble radius. 

Using the same assumptions as were used for the damped harmonic 
oscillator, the equation for the bubble volume pulsations shown in 
Figure 2 can be written as 



V = V(t) - Vo. 

The second-order differential equation for the oscillating bubble is 
then written as 

d^v dv 

m|^+p|^ + sv = 0. (3) 

The solution of equation (3) is the same as that of equation (1) 
and, therefore, the angular resonance frequency is still defined by 
equation (2), 



cod = cOo = Vs/m. 
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V(t) 




V = volume pulsation 
Vq s equilibrium volume 
V(t) = Instantaneous volume 



Figure 2 

Pulsating Bubble 

Source; Ref. 5:20 
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Bubble stiffness is defined by the change in pressure within the bubble 
versus the change in bubble volume as the bubble undergoes oscillation 
(Ref. 16], that is: 

s = - dp/dv. 

For a bubble containing gas at pressure Po, the bubble stiffness 
becomes 



where y = Cp/Cy is the specific heat ratio and is equal to 1.4 for air. 
The unit for the total pressure po is the atmosphere, which equals 
1.0133 X 10® dynes/cm2. The inertial constant m for a spherical 
bubble of mean radius Ro is defined as 



where p is the density of seawater (1.026 g/cm3). Substituting equa- 
tions (4) and (5) into equation (2) gives [Ref. 16], 



s = 7Po/Vo. 



(4) 



m = p/4tcRo 



(5) 




Using 



Vo = |itR3, 
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the final result for the angular resonance frequency of a gas bubble in a 
liquid as a function of bubble radius is 



CDo = 




or 



COo = 



1 . 

Ro V P • 



( 6 ) 



The resonance frequency fo is related to the angular resonance 
frequency coq by the relationship 

fo = coo/2tc. 

Therefore, the final result (in hertz) is 



Equation (7), alone, can be used to calculate the range of reso- 
nance frequencies for bubbles in seawater. More accurate equations 
which account for surface tension, viscosity, and thermal effects have 
been derived, but the effects of these factors on the resonance fre- 
quency of bubbles typical of those found in ship wakes is very small. As 
an example, the equation for the resonance frequency containing a 
surface tension term a is [Ref. 3], 




(7) 




( 8 ) 
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By examining equation (8). it can be seen that, as bubbles of increasing 
size are considered, the surface tension term becomes negligible (the 
value for surface tension in seawater is approximately 70.0 dynes/cm). 
However, surface tension and viscosity terms will become important in 
the next section when calculating the pressure amplitudes of the 
sound energy re-radiated from bubbles under the influence of dual- 
frequency excitation. 

Equations (7) and (8) show the relationship between bubble 
resonance frequency and radius. The most abundant bubbles in the 
surface ship wake environment are assumed to have approximate radii 
var3dng from 10 to 170 micrometers (|im) and extend downward from 
the surface to a depth corresponding to approximately three 
atmospheres of pressure [Ref. 171. Appendix A is a listing of a Fortran 
program for calculating resonance frequencies using both equations 
(7) and (8). For the purpose of comparison, the results of these two 
equations are tabulated in Tables la, Ib, and Ic using radii from 10 to 
170 mm and pressures from one to three atmospheres. These 
assumed values of bubble radii and their corresponding resonance 
frequencies are of primary interest for the Dual Frequency Pump 
method of bubble density measurement in ship wakes. 

B. DUAL FREQUENCY BUBBLE RESPONSE 

The theory of dual frequency bubble response is more complex 
than that previously used to find the resonance frequency of a bubble 
of a given radius. Part of this complexity is due to the inclusion of 
shear viscosity and surface tension terms, which become necessary for 
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TABLE la 



RESONANCE FREQUENCY VS. RADIUS FOR 
BUBBLES IN SEAWATER AT 1.0 ATMOSPHERE 



RADIUS 

(cm) 


RESONANCE FREQUENCY 
(KHz) 


FREQUENCY WITH 
SURFACE TENSION 
(KHz) 


0.0010 


324.1 


340.8 


0.0020 


162.1 


166.3 


0.0040 


81.0 


. 82.1 


0.0060 


54.0 


54.5 


0.0080 


40.5 


40.8 


0.0090 


36.0 


36.2 


0.0100 


32.4 


32.6 


0.0120 


27.0 


27.1 


0.0150 


21.6 


21.7 


0.0170 


19.1 


19.1 
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TABLE Ib 



RESONANCE FREQUENCY VS. RADIUS FX)R 
BUBBLES IN SEAWATER AT 2.0 ATMOSPHERES 



RADIUS 

(cm) 


RESONANCE FREQUENCY 
(KHz) 


FREQUENCY WITH 
SURFACE TENSION 
(KHz) 


0.0010 


458.4 


470.3 


0.0020 


229.2 


232.2 


0.0040 


114.6 


115.4 


0.0060 


76.4 


76.7 


0.0080 


57.3 


57.5 


0.0090 


50.9 


51.1 


0.0100 


45.8 


46.0 


0.0120 


38.2 


38.3 


0.0150 


30.6 


30.6 


0.0170 


27.0 


27.0 
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TABLE Ic 



RESONANCE FREQUENCY VS. RADIUS FOR 
BUBBLES IN SEAWATER AT 3.0 ATMOSPHERES 



RADIUS 

(cm) 


RESONANCE FREQUENCY 
(KHz) 


FREQUENCY WITH 
SURFACE TENSION 
(KHz) 


0.0010 


561.4 


571.2 


0.0020 


280.7 


283.2 


0.0040 


140.4 


141.0 


0.0060 


93.6 


93.8 


0.0080 


70.2 


70.3 


0.0090 


62.4 


62.5 


0.0100 


56.1 


56.2 


0.0120 


46.6 


46.9- 


0.0150 


37.4 


37.5 


0.0170 


33.0 


33.1 
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extremely small bubbles. Another reason for the complexity is that the 
bubble must now be viewed as a system of forced volume pulsations in 
the steady state. 

Once again, a bubble model must be utilized and certain 
assumptions made. The bubble shown in Figure 3 is being radiated by 
two sound fields (note that the bubble size is not scaled to compare to 
sound field wavelength). The bubble will remain approximately 
spherical throughout its volume pulsation and is considered to be sur- 
rounded by an infinitely extended incompressible fluid. Thermal 
effects and gas diffusion through the bubble wall are considered negli- 
gible. Also, damping of the bubble oscillation is considered through 
the inclusion of surface tension and shear viscosity in the bubble 
model with shear viscosity being taken as a constant for Newtonian 
fluids such as seawater. [Ref. 18] 

The mathematical form of the bubble model is a non-linear differ- 
ential equation for volume pulsations [Ref. 3], 

pRR + 3/^ pR2 = (po + 2a/Ro) (Ro/R)3y - (po + 2a/R) - 4pR/R + pi(t) 
+ Pp(t), (9) 

where 

Po = ambient pressure 

p = shear viscosity of the liquid 

a = surface tension 

Pj(t) = instantaneous pressure of imaging sound wave 
0)i = frequency of imaging sound 

Pp(t) = instantaneous pressure of pump sound wave 
(Op = frequency of pump sound 
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Liquid f ^ 
Properties \ p 




Imaging Sound 
Waves 




Figure 3 

Pulsating Bubble Undergoing 
Radiation by Two Soimd Waves 

Source: Ref. 18:283 
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Y = ratio of specific heats 

R = instantaneous bubble radius 

Ro = average bubble radius 

and the dots represent time derivatives. The incident sound waves 
are assumed to be'sinusoidal, that is, 

Pi(t)-= pi cos (Oit 
and 

Pp(t) = Pp cos (Opt, 

where pi and pp are simplified expressions for the peak pressure 
amplitudes (imaging and pump, respectively) at a particular distance 
from the sound source. 

In order to solve equation (9), it is useful to expand the instanta- 
neous radius R in a Taylor series, 

R = Ro (1 + x) = Ro + xRo, (10) 

where Ro is the mean bubble radius and | x | is much less than one. It 
is emphasized that equation (10) is valid only for small volume pulsa- 
tions in the linear region. Because the pulsations are small, and 
because the solution for the dual-frequency excitation method 
requires only first and second harmonics, cubic and higher order 
terms of the expansion are ignored [Ref. 19]. Substitution of equation 
(10) into equation (9) yields. 
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pR2(l + x)x + 3/2 pR^x2 = (po .+ 2a/Ro)[l - 3tx + 3/2 7(87 + l)x2] 

- [po + (2a/Ro)(l - X + x2)] - 4|x(x - xx) + pi cos coit 

+ Pp cos COpt. (11) 

Equation (11) is now a second -order differential equation for volume 
pulsation in terms of x. 

The Dual Frequency Pump Method involves incident sound waves 
at two different frequencies. Therefore, a solution for x using dual 
frequency excitation includes the first harmonics (coi and cop), the 
second harmonics (2 cdi and 2cop), and the two sidebands (coi + cop and 
coi -cop). With all frequencies included, the solution may be written as, 

X = Ao + Ai cos (coit + 0i) + A2 cos (copt + 02) 

+ A3 cos (2o)it + 03) + A4 cos (2copt + 04) 

+ A5 cos [(coi + cop)t + 0s) + Ae co^ [(coi - cop)t + 0e). ( 1 2) 

The solution for x contains the amplitudes Aj for the changes in radius 

during pulsations. For example: 

Ao = amplitude change of the average radius, 

Ai = first harmonic amplitude change for coi, 

A2 = first harmonic amplitude change for cop, 

A3 = second harmonic amplitude change for coi, 

and so on. The phase angles 0j are the relative phase differences 

between the incident wave and the particular harmonic pulsations. 

[Ref. 3, 19] 
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The goal of this development is to derive a set of equations giving 
the pressure amplitude in seawater as a function of the distance r from 
the bubbles insonlfied by the Dual Frequency Pump Method. These 
pressure amplitudes, when measured at discrete frequencies, can be 
used to determine the number of bubbles at each radius. The tools to 
do this have now been supplied. The amplitudes of the radius excur- 
sion, Aj, may be solved for by substituting equation (12) into equation 
(11). Algebra for this step is quite lengthy, so only the results are 
shown below. The problem-solving technique is simpler for a single 
sound wave (Miller’s single sound wave reduces the assumed solution 
from seven terms to three) [Ref. 19]. Numerical values for the ampli- 
tudes of the radial pulsations in equation (12) are best solved with a 
computer using the following: 



Ai = PiXi/paX 




(13) 


A2 = PpX2/pO)|R2 


- 


(14) 


A5 = A1A2X22X5 




(15) 


A6 = AiA2X’i2X6 




(16) 



where, 

(17) 



Xi = [(1 -nf)2 + 52 Q 2 J- 1/2 

X2 = [(1 - n|)2 + 52Q2J-1/2 



(18) 




3 _ ^ 3y (3y + l)(po + 2a/Ro) - 4a/Ro^ 
+ ^ niQ 2 + 2pR^co^ 



2 



( 19 ) 




3 ^ _ 3y (3y + l)(po + 2a/Ro) - 4a/Ro^ 

- 2^1 ^2 + 2 p ^ 



2 



( 20 ) 



X5= [d -Qi)2 + 52 Q 2 ]-l /2 (21) 

X 6 = [(1 - Q |)2 + 52q2 ]-1/2 . . ( 22 ) 

The parameter Q is a dimensionless ratio relating the pump, image, 
and sideband frequencies to the resonance frequency coq of the bubble. 
That is. 



Ql = COi/COo = COp/cOo 

Q 5 = (o)i + cop)/a)o Cle = (coi - cop)/o)o (23) 

The expression 5 is the viscous damping coefficient without thermal 
cind radiation considerations. It is written as 

5 = 4p/pcooRj. 

The radiated pressures (fundamental and sideband) from the bubble 
excited by two frequencies may now be determined in terms of Aj 
using equations (13) through (23) [Ref. 3]. 
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In order to derive the equation for pressures as a function of the 
distance from the bubble, it is necessary to model the bubble as a 
monopole radiator [Ref. 14]. A monopole radiator is really nothing 
more than a simple source. The equation for pressure at a distance r 
from a simple source is [Ref. 14] 



Since the Dual Frequency Pump Method requires that only amplitude 
be considered, and not phase, the above equation can be rewritten as 



Equation (24) represents the pressure amplitude for the frequency. 
Recall that 




Pj (r,t) = Pj eJ(“t - 



(24) 



where the pressure amplitude is 




(25) 




(26) 



and for the pulsating sphere [Ref. 14] 



Q = 4tcR2Uo 



(27) 
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where Uq is the speed amplitude. In order to solve equation (25) using 
amplitude changes, equation (27) must be in terms of the mean radius 
Ro and the amplitude changes It is shown in Appendix B that 

Uo = RoAjCOj. 

Therefore, 

D _ 

4%r 

47cR2UoC0j 
“ Po 47ir 

which 5 delds 



Pj = poRjAjCoJ/r. (28) 

Equation (28) is the pressure amplitude at a distance r from a bubble 
under the influence of a sound wave at frequency coj. For the Dual Fre- 
quency Pump, the frequencies are coj = coi, cop, and coi ± cop. The 
amplitudes are Aj = Ai, A 2 , As, and Aq. Equations (13) through (22) 
are substituted into equation (28) to give the desired equation for 
pressure amplitude that results from the imaging, pump, and sideband 
energy reradiated by the bubble. These pressure amplitudes are Pi, 
P 2 , P 5 , and Pe, respectively, where P 5 and Pe will be referred to as the 
upper and lower sideband pressure amplitudes P+ and P_, respectively. 
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-The key equations are summarized below for Pi, P2, P+, and P_: 



Pi = Pi Xi nf Ro/r 



P2 = pp X2 

P+ = X12 X5 (coi + i0p)2 Rj/r 



P- = X'i3 Xe (coi - cop)2 R3/r. 



V, L. Newhouse and P. M. Shankar simplified these equations because 
the imaging frequency is much greater than the pump frequency. It is 
possible to say 



X5 = Xi 



X6«Xi 



and 



X12 = [Ref. 3] 

The advantage of dual frequency bubble detection can be demon- 
strated by examining the sum sideband pressure amplitude P+. 
Consider a single bubble of resonance frequency coq. When the pump 
frequency sweeps past the bubble’s resonance frequency (cop = coq), the 
bubble undergoes maximum oscillations. At this point, the sum fre- 
quency sideband pressure is at its peak and can be written as 
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(28a) 



o PPiPp Ro 1 
~ (pcooRo)-^; S ■ r • 

At all other frequencies, that is cop coq. the expression for P+ becomes 



P+ 



_ PPlPp Y ^ 
” (pcooRo)^ ’ ^ ‘ r * 



(28b) 



It should be noted that the X 2 term of equation (28b) contains the 5 
term of equation (28a) and the dimensionless frequency term Q2- 
When cop = (Dq, the CI 2 term is equal to one and equation (28b) 
becomes equation (28a). In equation (28a), with the pump frequency 
equaling the bubble resonance frequency, the pressure amplitude of 
the sideband is proportional to the mean bubble radius Rq. At all other 
frequencies, the upper sideband pressure amplitude is proportional to 
both X 2 and Rq as shown in equation (28b) . This is significant because 
the sideband pressure P+ radiated by the bubble provides a maximum 
pressure amplitude for a particular bubble radius, even though that 
sideband pressure amplitude is much smaller than that of the pump 
pressure amplitude P 2 . For a single bubble. Figure 4 shows the 
correlation between the sideband pressure amplitude P+ of curve “a” 
and the pump pressure amplitude P 2 of curve “b” when plotted as a 
function of the pump frequency fp in kHz. Figure 4 demonstrates that, 
for a single bubble, the maxima, or spikes, of both pressure amplitudes 
coincide with the resonance frequency of the bubble. (Ref. 3] 

If the sideband pressure amplitude P+ and the pump pressure 
amplitude P 2 are plotted as a function of the mean bubble radius, the 
difference between single-frequency resonance excitation and the 
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difference between single-frequency resonance excitation and the 
dual-frequency pump method becomes clear. Figure 5 shows the 
pump pressure amplitude on curve “b” and the sideband pressure 
amplitude on curve “a” plotted against the mean bubble radius. This 
means Figure 5 considers bubbles of many radii with the frequency of 
the pump sound field fixed. While the pump pressure amplitude has a 
local maximum pressure at the point where bubbles of corresponding 
resonance frequency (cop = coq) are emitting maximum energy, the 
pump pressure amplitude rises above that local maximum as energy is 
reflected from the larger, nonresonating bubbles. The sideband pres- 
sure amplitude P+ of curve “a,” however, still maintains a single sharp 
peak for the one bubble radius undergoing resonance and does not 
show any sound energy refiected from larger bubbles. Sound pressure 
from the nonresonant bubbles is not present in the sideband pressure 
peak. The reason for this difference lies in the fact that sideband 
production is a nonlinear process which occurs only when • the 
pulsation amplitude is large. In general, only resonant bubbles 
undergo large enough pulsations to produce the sidebands. Larger, 
nonresonant bubbles scatter significant amounts of the pump and 
imaging fields, but their pulsation amplitude is too small to produce 
the sidebands. Thus, the Dual Frequency Pump method can distin- 
guish bubbles of different sizes much better than single -frequency 
excitation. [Ref. 3] 
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Figure 4 

Pressure Amplitude (P^ and ) vs. Frequency (ly 
Source: Ref. 3:1475 




XIO^ 



0.3 



0.2 




0.1 



Figure 5 

Pressure Amplitude (P^and PJ vs. Radius (Rq) 



Source: Ref. 3:1475 
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Wien bubbles of different radii are irradiated with both the pump 
and imaging fields, several “sum-difference” sideband pressure peaks 
corresponding to those various bubble radii will be present. The mag- 
nitudes of these pressure peaks give an indication of how many 
bubbles of a particular radius are present. The methods of determin- 
ing the number of bubbles of a particular radius will be discussed 
further in the “Considerations” section of the next chapter. 

C. TRANSDUCER BEAM PATTERNS 

It is now necessary to consider the volume of water insonified 
while sampling the bubble cloud. This sample volume is important in. 
determining the bubble cloud density. Because the sample volume is 
specifically determined by the beam patterns of the imaging and pump 
transducers, the final theory section will study far field beam patterns 
and their importance. 

The first step in defining the sample volume is to examine the 

near and far field of the transducer’s beam pattern. Sound waves in 

I 

the near field have pressure amplitude “nulls” and phase shifts as a 
function of the range which are undesirable when using sound to make 
measurements. Sound waves in the far field are expanding spherically 
and the pressure amplitudes decrease only as a function of range from 
the sound source. There are no pressure amplitude “nulls” or phase 
shifts in the far field. Therefore, it is best to have the target of inter- 
est in the far field, where calculation of the sound pressure amplitudes 
is easiest. Several methods of calculating the start of the far field are 
available. Here, a conservative estimate of the far field range will be 
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lerived. Also, the angular beam widths are needed to define the 
5 ample volume. These angular beam widths and far field points for 
Doth pump and imaging sound field transducers wiU be calculated for 
representative frequencies used by the Dual Frequency Pump method. 

The Dual Frequency Method of bubble detection is based on 
sensing the radiated pressure amplitudes of the bubbles at the side- 
Dand frequencies of the imaging sound field (coi ± cop). In calculating 
die transmitted pressure amplitudes at the point where the bubbles 
are located, both the pump and imaging sound waves are considered 
^0 be planar at that location. Insonilying the target bubbles with sound 
ivaves independent of phase and without nulls requires that the 
bubbles be in the far field of the transducer beam. To find expressions 
for the start of the far field for both the pump and the imaging trans- 
ducers, it is easiest to take a linear systems approach. The physical 
situation can be described as a volume-type aperture serving as one of 
the transducers shown in Figure 6. A volume aperture is a description 
which can apply to any shape of sound source which is transmitting 
acoustic signals. Figure 6 shows the vector r from the coordinate 
origin to the target point in space. The vector r© is from the coordi- 
nate origin to a point on the surface of the volume aperture and iden- 
tifies the location of each sound source. The volume aperture can be 
either the pump or imaging transducer which transmits sound energy 
or, due to reciprocity, can be the receive transducer which senses the 
radiated sound energy from the bubbles. The velocity potential 
solution to the wave equation can be written in terms of the free space 
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Green’s function. It is from Green’s Function that the near and far 
fields of the transducer beam pattern are defined. [Ref. 20] 

The free space Green’s Function can be written as 



g(r/ro) = 



exp [-jk|r -rj] 



r -ro 



where 



k = 2 tc/X. 

Here, the I r - r o I range term appears as both an amplitude 
(denominator) and a phase (exponential). To determine the point 
where phase is no longer a factor, it is necessaiy to define the maxi- 
mum extent of the near field. Green’s Function is used for this pur- 
pose, and after undergoing a binomial expansion can be written as 

g(r/ro) « H^r) ^ [jk(ar • Tq)] exp [-Jk(r2 - (ar • r.o)2/2r)] (29) 

where ar is a unit vector in the direction of r . The last complex term 
of equation. (29) is the critical expression in determining the near 
field. It is written as 



exp [-Jk(r2 - (ar • ro)2/2r)] . (30) 

The magnitude of the exponential argument in equation (30) is 
written as 



% [rJ - (ar • Lo)2]/Xr 
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Figure 6 

Volume Apertiure Transducer 

Source: Ref. 20:2 
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and can be rearranged to yield 



7i[r2 - (ar • Tq)^]/X 



( 31 ) 



r 



Equation (30) will be a significant phase term of the free space 
Green’s Function if equation (31) becomes large. Equation (31) 
becomes large if the numerator is larger than the range r. Because it 
is defined as the region having significant phase, equation (31) can be 
used to define the near field. Phase is significant when 



A Fresnel approximation is now used, which is equivalent to setting 
the dot product in equation (32) equal to zero, that is [Ref. 20] 



If the worst case is considered, the radial aperture distance r© can be 
replaced by the maximum radial dimension of the aperture R. The 
final form of equation (32), the extent of the near field, is now written 
as 



r < 7t 



[r| - (Sr ■ to)2)/X. 



(32) 



(^r ■ to)^ 0- 



r < 7 cR2/X, 



(33) 



where 



r = rainge from the transducer 
R = maximum radial dimension of the transducer 
X = wavelength of transmit signal. 
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Further manipulation of the transmitted acoustic signal with Fourier 
Transforms )delds the Near Field Directivity Function D(f,r, cd in terms 
of the complex aperture function (At), spatial frequencies fet), and 
range (r) (Ref. 20]. It is 



The Near Field Directivity Function describes a beam pattern which 
does not assume plane waves over a small space since the directivity is 
a function of the range and phase terms. 

As previously mentioned, it is desirable to be in the far field, 
where the directivity is not a function of range. The far field point is 
now defined as 



The Far Field Directivity Function is obtained by using Green’s Func- 
tion and a Fraunhofer approximation [Ref. 20]. The result is 



Equation (36) is not a function of range and therefore will be used 
with equation (35) to develop transducer placement for the Dual Fre- 
quency Method of bubble detection. While it is not important to 
understand all the terms of equations (34) and (36), it is important to 
see that the acoustic frequency, the physical size of the transducer. 



Dt (f,rj^ = j At (fio) e-Jk^o/2r ^2 k (^- r^) dVo 



(34) 



r > 7 iR2/X. 



(35) 




(36) 
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and the range to the target point from the transducer are key factors 
in determining beam patterns. 

Previously, the general physical situation was that of a transducer 
considered to be a volume aperture (Figure 6) from which the general 
form of the Far Field Directivity Function (equation (36)) was obtained. 
It is now possible to be more specific. The transducers to be used in 
the Dual Frequency Pump method of bubble detection are closely 
approximated by planar, circular piston apertures. The circular piston 
aperture is shown in Figure 7. The general expression for complex 
aperture function is written in polar coordinates as 

A(f,r,0) = a(f,r,0) . 

Circular symmetry is assumed for both the above aperture function and 
the Far Field Directivity [Ref. 20]. Therefore, they can be written as 



D(f,0,\|/) = Ho (Ar (f,r)), 

where Ho is the zero order Hankel Transform [Ref. 20]. The Directiv- 
ity Function is independent of the angle due to circular symmetry. 
Because the Hankel Transform can be written as 




and 



Ho [Ar (f,r)) = 2k 



f rAr(f,r)J 




o 
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Figure 7 

Circular Piston Aperture 

Source: Ref. 20:80 



53 



the Far Field Directivity Function of the circular piston is 



D(f,0) = 2k f rJo l^^^sln 0^ dr 



where “a” is the radius of the transducer face as shown in Figure 7. 
Using the identity for Bessel Functions 



aJo (a) da = xJi[x], 
o 



the Directivity Function is now written as 



D(f,0) 



J 1 \[2k3./X) sin 0] 
^ sin 0/X 



(37) 



where 

a = radius of the transducer 
0 = angle off of transducer axis (half beamwidth) 

X = wavelength of transmitted sound 
Ji(x) = first-order Bessel function. 

The normalizing factor for equation (37) is given by the directivity 
along the transducer axis, and can be written as 



D(f.0) I e=o = 7ca2. 



Therefore, the Normalized Far Field Directivity Function is 



_ 2Ji[(27ca/X) sin 0] 

DnU,0) - (2K2i/X) sin 0 



(38) 
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In order to define the volume of seawater that will be insonified by 
the transducers, equation (38) must be used to determine the 
beamwidths of the main lobes at their 3 decibel down points. No 
standard value for the ratio of intensities has ever been established for 
determining the beamwidth, but 3 decibels has been used by many 
authors and will be used here [Ref. 14]. Letting 



20 Logic Dn (f.e) = -3 dB (39) 

so that 



Dn (f,e) = = .7079 



then equation (38) becomes 



2 Ji [X] 

X 



= .7079 



by substituting 

X = [(2KQ./X) sin 0]. (40) 

A solution can be found for x by referring to a standard table of 
first-order Bessel Functions for a circular piston (Ref. 14). The result 
is X = 1.614. The half angle 0 of the beams can now be calculated for 
each transducer using equation (40). The results of the far field 
determinations (equation (35)) and the beamwidths are given in Table 
II. The sound speed in seawater was assumed to be 1,500 
meters/second. This information provided a basis for determining 
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transducer placement and calculating the volume of seawater which 
contained the target bubbles. A total of three transducers were used 
during the experiment. Two high-frequency transducers provided the 
imaging sound field and the receiver, and the third transducer pro- 
vided the pump sound field. Table II results are for pump frequencies 
of 100 and 500 kHz, and an imaging frequency of 2.25 MHz to reflect 
samples of the frequencies used. 
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TRANSDUCER BEAM PATTERN DATA 




57 



in. PROBLEM APPROACH 



A. OVERVIEW 

As discussed in the introduction, the goal of this thesis experi- 
ment is to show that the Double Frequency Pump technique is practi- 
cal for ship wake bubble density measurements. Chapter II provided 
the basic theoretical background for calculating the resonance fre- 
quency of different size bubbles, the pressure amplitudes at the sum- 
difference sidebands, and the far-field beam patterns of the 
transducers. 

Other factors to be considered in the practicality of this bubble 
detection technique include the frequency range of the swept pump 
sound field, the sideband pressure amplitudes, sample volume place- 
ment and sizes, distinctive interference, bubble screening, transducer 
orientation, and transducer response. The statistical sample time that 
a point in the bubble cloud is measured and the imaging frequency are 
also important. 

B. CONSIDERATIONS 

1. Pump Frequencies 

Knowing the range of resonance frequencies for the different 
bubble sizes present is essential so that the correct pump frequencies 
can be used. In order for the bubbles to radiate the maximum amount 
of sound energy at the dual-frequency sidebands (coi ± cOp), the bubbles 
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must be driven at resonance. The bubble sizes chosen here are those 
ranging in radius from approximately 10 to 170 micrometers (|im). 
Bubbles much smaller than 5 |im in radius resonate at very high fre- 
quencies and do not affect sound of particular interest to the surface 
ship wake problem at the ocean surface. Due to the effects of diffusion, 
small bubbles diffuse and disappear into the. water. Larger bubbles (above 
the 170-200 |xm radius range) tend to rise quickly or split up into 
smaller bubbles, and therefore do not persist very long in the wake. 

In order to drive the bubbles present in the bubble cloud at 
resonance, the selection of the pump frequency range and the 
pressure amplitude (pp) is important. However, these characteristics 
also control the harmonics produced as the bubbles react to the pump 
sound field. Harmonics of the pump sound field and the resonating 
bubbles are critical because they can create background noise at 
higher frequencies. If the pump field frequencies (cop) are of a very 
wide range, the harmonics (each of wide frequency) extend farther up 
the frequency spectrum. If the pump field amplitude is too high, the 
harmonics created by bubble resonance will be multiple, and again the 
result is the same. Harmonics can defeat the detection of the dual- 
frequency sidebands about the imaging frequency if the resonance 
harmonics extend to high enough frequencies. If the pump frequency 
range is small, the harmonics induced will have a short frequency 
range and, therefore, will not extend to higher spectrum frequencies. 
If the pump pressure amplitudes (pp) at the bubble location are kept 
low (yet large enough to resonate bubbles), the number harmonics will 
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be less and the result the same; the harmonics will not extend to 
higher frequencies on the spectrum. 

2. Sideband Pressure Amplitude 

The pressure amplitude of the sum sideband (P+) that is cal- 
culated using equations (28a) and (28b) makes it possible to deter- 
mine the exact number of bubbles of a specific radius present within a 
given sample, or measurement, volume. This is done by calculating 
the radiated pressure amplitude (P+) for each bubble radius present 
and then comparing the results with the actual received pressure 
amplitudes. Table III contains sample results from equations (28a) 
and (28b). Appendix C contains a listing of a simple Fortran program 
which will calculate the P+ amplitudes for bubbles using equations 
(28a) and (28b). The Fortran program uses input values for pump 
pressure amplitude (pp) and imaging pressure amplitudes (pi) that are 
calculated in Appendices D and E. The bubbles were assumed to be at 
a distance (r) of 6.7 cm from the transducer. An example of the output 
at one atmosphere of pressure is plotted in Figure 8 for a single 
bubble. The curve that has been plotted is not actual data, but only 
represents the pressure amplitude as it rises from the noise. The data 
points are the computer results. Figure 9 is plotted for bubbles of the 
same radii as that shown in Tables I and III (same curve estimates as 
Figure 8). For the sideband pressure calculation to be accurate, the 
location of the target bubbles relative to the transducers must be 
known. Also, the size of the volume in which the bubbles are located 
is important. This volume of water where the bubbles are detected 
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TABLE III 



THEORETICAL CALCULATION RESULTS FOR 
UPPER SIDEBAND PRESSURE AMPLITUDES 



BUBBLE 

RADIUS 

(cm) 


RESONANCE 

FREQUENCY 

(kHz) 


SUM SIDEBAND 
PRESSIJRE PEAK 
(dynes/cm2) 


SUM SIDEBAND 
PRESSURE NOISE 
(dynes/ cm2) 


0.0010 


324.1 


0.0010 


0.000429 


0.0020 


162.1 


0.0040 


0.000925 


0.0040 


81.0 


0.0161 


0.001890 


0.0060 


54.0 


0.0362 


0.002846 


0.0080 


40.5 


0.0643 


0.003801 


0.0090 


36.0 


0.0814 


0.004277 


0.0100 


32.4 


0.1004 


0.004754 


0.0120 


27.0 


0.1446 


0.005707 


0.0150 


21.6 


0.2260 


0.007135 


0.0170 


19.0 


0.2903 


0.008088 



61 




OOONCDlCTtCOcNrHQ 

S S B S S S S S S S 

(^3/S9tLCp) 9pn;ndmv 9Jnss9Jd pirBq9pis 



62 



Pump Frequency (KHz) 

Figure 8. Plot of Computed Sideband Sound Pressure Amplitude 
For a Single Bubble of Sixty Micron Radius 
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Bubble Radius (cm) 

Figure 9. Computed Sideband Sound Pressure Amplitudes 
For Bubbles of Many Different Radii 



and measured will be called the sample volume. Sample volume con- 
siderations are discussed in the next three sections.. 

3. Sample Volume Placement 

The first parameter which could define a sample volume is 
the Far Field of the beam pattern. As mentioned in the section of 
Chapter II on Transducer Beam Patterns, it is considered desirable 
that the target bubbles lie in the far field of the beam pattern to avoid 
the “nulls’ and phase shifts found in the near field. In the far field, 
the radiated pressure amplitude of a sound source drops steadily pro- 
portional to the inverse of the range [Ref. 14]. The equations for cal- 
culating the sum-difference pressure amplitudes (28a and 28b) 
require that the pump and image sound pressure amplitudes be 
known at the point where the target bubble is located. It is easiest to 
calculate the pump and Imaging sound energy present at the location 
of the bubbles in the far field. Due to the principle of reciprocity, 
which now considers the target bubble as the sound source, it is also 
desirable that the target bubble be in the far field of the receive 
transducer so that the sideband pressure amplitude P+ may be accu- 
rately calculated at the receive transducer face. 

Placement of the target bubbles in the far field of the beam 
pattern appears to be ideal. However, it is important to provide 
enough imaging and pump sound pressure to the bubbles to create the 
sideband pressure amplitudes and raise them above the noise level. 
Dr. P. M. Shankar noted that it may be necessary to locate the imaging 
and pump transducers close to the bubble cloud, that is in the near 
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field, to provide the bubbles with these pressures [Ref. 21]. The 
receiver may also need to be close to the bubble cloud to detect the 
sound pressures reflected and radiated by the bubbles. It is possible to 
make this argument in favor of a near field sample volume placement, 
despite the fact that it complicates the calculations of the imaging and 
pump pressure amplitudes (pi and p 2 ) presented to the bubbles. The 
effects of transducer placement closer to the target bubbles will, 
therefore, be examined. 

4. Destructive Interference 

Consider the example of two bubbles of exactly the same 
radius within a sample volume where all transducer beams (Pump, 
Image, and Receive) meet. It is possible for destructive interference 
to occur due to differences in the range from the receive transducer 
while the bubbles are emitting sound energy at the sideband 
frequency, coi + cop. This results in the possibility of not recording the 
full P+ amplitude for two bubbles in the sample volume at the same 
time. The likelihood of destructive interference can be reduced by 
designing a small or a thin sample volume in order to reduce the 
probability of two or more bubbles of the exact size being present in 
the sample volume at the same moment. 

5. Bubble Screening 

Bubble screening could also cause a false recording of the 
sideband pressure amplitude for a particular bubble. Bubble screening 
occurs when a large bubble blocks the radiated energy from a smaller, 
resonating bubble [Ref. 12]. Again, a possible solution would be the 
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design of a small or t±iin sample volume so that the receive transducer 
has a clear view of all bubbles passing through the sample volume. 
Fewer bubbles in the smaller volume cause less screening. Having 
bubbles present outside the sample volume would also screen the 
bubbles as they radiated sideband sound pressures. 

6. Sample Volume Size 

The factors of far field, destructive interference, and bubble 
screening all favor the use of a small sample volume to measure bubble 
density at points within a bubble cloud. These factors also favor 
isolating the sample volume from the bubble cloud. The bubble cloud 
could screen the sound from the sample volume. The placement of 
the imaging and receiving transducers can be used to determine -the 
sample volume size. Table II provided values for the far field points 
and beam widths for the 2.25 MHz imaging and receive transducers. 
Using the Table II values produces a calculated sample volume of 
approximately .107 cm3. The calculations and assumptions are shown 
in Appendix F. This sample volume may to be small enough to 
eliminate substantial bubble screening or destructive interference. 

7. Transducer Orientation 

Another important consideration which determines trans- 
ducer placement is the direction that the transducers are pointed. It 
is undesirable to directly radiate; the face of the receive transducer 
with either pump or imaging sound fields. Direct irradiation increases 
the probability of raising the noise level and obscuring the desired 
dual-frequency sidebands. In addition, direct irradiation can result in 
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the generation of sum-difference sidebands. In the Shankar and 
Newhouse experiment, the image and receive transducer beams 
intersected at right angles. The pump transducer was set slightly off 
the vertical in a downward-looking configuration. The pump trans- 
ducer orientation used by Shankar and Newhouse was compared to an 
upward-looking transducer in our experiment. 

8. Transducer Frequency Response and Side Lobe Effects 

Another factor in considering the sound pressures provided 
to the bubble by the imaging and pump transducers is the frequency 
response of each transducer. The section on transducer beam 
patterns in Chapter II pointed out that the output of the transducer is 
dependent on the input signal. If the frequency of the input signal is 
varied, the acoustic output at a fixed distance from the transducer face 
will also vary depending on the frequency response of that particular 
transducer. This is not a problem for the imaging transducer as the 
image frequency coi is fixed. The acoustic output is, therefore, also 
fixed. However, it is a problem for the pump transducer as the pump 
frequency cop is swept. The pump transducer will not provide flat 
sound pressure amplitude values at a fixed point in space when the 
pump frequency is swept and, therefore, does not provide a fixed value 
of pp to be used in equations (28a) and (28b). Some examples of the 
different sound pressure levels for the pump transducer frequency 
response and calculations of the pump pressure amplitude (pp) from 
those levels are provided in Appendix D. Shankar and Newhouse also 
encountered this frequency response problem in their original Dual 
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Frequency bubble detection work [Ref. 12]. They solved the problem 
by slowly stepping through the pump frequencies and adjusting the 
power of the input electrical signal to the pump transducer. Adjusting 
the input electrical power at the different frequencies gave a fairly flat 
frequency response of the pump sound field and allowed the 
calculations of equation (28a) and (28b) to be carried out [Ref. 21]. 
This method of obtaining flat frequency response is fine in the 
laboratory, but would prove difficult due to the bubble detection 
problem in surface ship wakes. Surface ship wakes generate bubble 
clouds of unknown size distribution; therefore a constant acoustic 
pressure level at all frequencies is necessary. One of four methods 
might be available to obtain the correct output. Tlie first might be to 
control the input level to the pump transducer with a computer so 
that as the frequency range of interest is swept the power is adjusted. 
A graphic equalizer performs the same function. The second might be 
to find a high-quality transducer that has a flat frequency response 
over the range of pump frequencies used. The use of several pump 
transducers of different center frequencies to cover the pump 
frequency range may also work. Third, a high-quality, calibrated 
broadband noise source may provide equal levels of acoustic energy to 
the sample volume at all the pump frequencies. Lastly, a computer 
could solve the problem by calculating the sideband pressure 
amplitudes for every different input pump pressure amplitude. 

Yet another source of inaccurate bubble measurement using 
acoustic means could be the energy transmitted or received by the 
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side lobes of the transducers. It is not likely that side lobes are a 
problem because the first side lobe of a circular planar transducer is 
approximately 17 decibels down from the main lobe [Ref. 20]. Due to 
the relatively small amounts of energy available from the resonating 
bubbles, side lobes are probably not a factor in the accuracy of bubble 
detection or the transducer placement. 

9. Statistical Sample Times 

Assuming that the Dual Frequency Pump Method detects and 
counts the bubbles passing through the sample volume, it should also 
be considered as to how many sample volumes must be “looked at“ by 
the receiver to consider the bubble density measurement accurate. 
The number of sample volumes “looked at” can be converted to the 
time duration of the “look.” if the speed of either the measurement 
device or the bubbles passing through the device is known. 

The number of sample volumes necessary to make the bubble 
density measurement statistically accurate for a given bubble stream is 
given by a general sample size equation [Ref. 22]. This equation is 
written 






^/2 
E J 



(41) 



where. 



n = number of samples 

Za /2 = confidence factor from the normal standard distribution tables 
E = maximum allowed error 



69 



Once the number of sample volumes necessary is determined by using 
equation (41), the duration of the sample is calculated using the speed 
of the bubbles and the dimensions of the sample volume, specifically 
the sample volume length. Appendix G provides an example of this 
calculation technique using a bubble rise speed of 1.50 cm/s [Ref. 23]. 
The statistical sample duration, for example, needed to provide a 95- 
percent confidence factor and only 5-percent error is approximately 
132 seconds. 

C. EXPERIMENTAL PROCEDURE 

This study of dual frequency bubble detection was carried out in 
three phases. The first phase was conducted in a small ten-gallon 
aquarium. It used fresh water for the purpose of making initial 
measurements, adjustments, and calibration. The second phase was 
also conducted in the ten-gallon aquarium, but this time using 
seawater. This phase was used to test the seawater medium and to 
test or refine different transducer arrangements. The third and final 
phase was carried out in a large, acoustically insulated fresh-water 
tank to investigate the effects of an acoustically quiet environment on 
bubble detection. The large insulated tank measured approximately 
four by two by three (4x2x3) meters. 

As mentioned previously, the first two phases of the experiment 
for the Dual Frequency Pump Method were done in a ten-gallon aquar- 
ium. This experimental setup is shown in Figure 10. The aquarium was 
first filled with purified deionized water produced in the laboratory, 
and later with seawater drawn from the Hopkins Marine Station in 
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Figure 10. Experimental Aquarium Setup and Equipment Rack 







Monterey. The seawater was relatively clean and closely approximated 
the salinity of open ocean seawater. 

Bubble generation techniques were briefly investigated with elec- 
trolysis being selected as the best method. Electrolysis provided for 
easy and controllable bubble clouds and generated bubbles of the size 
closest to those found in surface ship wakes. Insulated brass wire with 
the ends exposed was found to produce the fewest by-products and 
the best bubbles. The brass wire was laid on the bottom of the 
aquarium and the bubbles were allowed to rise through the sample 
volume of the transducer beam patterns under their own buoyancy. 
Power for bubble generation was supplied by an HP 6237B Power 
Supply and varied between .5 and 10 volts. -Both large bubble clouds 
and small bubble streams were generated by varying the voltage to test 
the effects of bubble screening. 

The transducers were fixed within the aquarium (AQ), and later in 
the large acoustic tank (LT), using plexiglass mounts. Three transducers 
were used during the entire experiment. For testing in the aquarium, 
both the imaging (center frequency 2.25 MHz) and receive (center fre- 
quency 2.25 MHz) transducers were mounted in the same plane at a 90- 
degree angle to each other as shown in Figures 10 and 1 1. This transducer 
mount had two sets of holes available for mounting the image and receive 
transducers such that the sample volume (defined by the beam inter- 
section) could lie in either the near (image near field or INF) or the far 
fields (image far field or IFF) of the beam patterns. The pump transducer 
(500 kHz center frequency), which was physically the largest of the 
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Figure 11. Image, Receive, and Upward Looking Pump Transducer Mounts 




t±iree, was mounted with two different placements so as to study the 
effects of radiating the rising bubble cloud at an upward or downward 
angle. The upward-facing pump transducer mount (designated “UP” 
for Upward Pump) for use in the aquarium is also shown in Figures 1 1 , 
12a, and 12b. The downward-facing pump transducer arrangement 
(designated “DP” for Downward Pump) can be seen in Figures 13, 14, 
and 15. For testing within the aquarium, the pump transducer mounts 
were kept separate to allow flexibility in movement of the mounts 
while investigating the best transducer alignment relative to the bub- 
ble stream. Testing during phase three (in the large tank) required that 
the transducer mounts be a single apparatus as shown in Figure 15. 

AH the transducers used in our experiments were manufactured 
by Panametrics Inc. The high-frequency imaging and receive trans- 
ducers had a center frequency of 2.25 MHz. The low-frequency pump 
transducer had a center frequency of 500 kHz. The output of all three 
transducers was checked using a calibrated hydrophone. 

The receive transducer was used to pick up signals reflected or 
radiated from the bubbles, or to receive noise from the other trans- 
ducers. The signal received by the receive transducer was amplified 
by a 27 dB pre-amplifier. The signal was then passed on to an 
HP3585A spectrum anal)^er and a Kikusui C056100A oscilloscope. 
The oscilloscope was used for several applications, including the 
measurement of reflected signals, calibration checks of the signal 
generators, and as an indicator of the generated bubble stream 
placement. The HP3585A spectrum analyzer was used for picking out 
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Figure 12a. Imaging (Far Field), Receive, and Upward Looking Pump (AQUPIFF) ih Aquarium 
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Figure 12b. Imaging (Near Field), Receive, and Upward Looking Pump (AQUPINF) in Aquarium 
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Pump Transducer Mounts 
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Figure 14. Imaging (Far Field). Receive, and Downward Looking Pump 
(AQDPIFF) ih Aquarium 
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the primary and sideband sound pressure levels in the frequency 
domain. The limiting equipment in the receive process was the 27 dB 
pre-amplifier, which had a 50 kHz low-frequency cut-off. The 50 kHz 
low-frequency cutoff was not a problem when measuring dual- 
freqeuncy sideband levels, but it was a complication when measuring 
low-frequency transducer responses or fundamental bubble resonance 
levels. 

The image transducer was driven by a HP 3314A Function 
Generator. This function generator could provide imaging signal fre- 
quencies up to 10 MHz and input levels up to 10 volts zero to peak 
with little distortion. 

The pump transducer was provided input signals from a Wavetek 
Model 22 Sweep Generator. The signal was amplified ten times (x 10) 
through a HP 467A power amplifier. The amplifier was added after a 
conversation with P. M. Shankar of Drexel University revealed the 
. need for sweeping the pump frequency range with higher pressure 
amplitudes than were originally used [Ref. 21]. All the equipment used 
during the experiment is pictured in Figure 10. Data from each phase 
was recorded and plotted from the spectrum anal3^er using a HP 7090 
plotter. 

1. Phase One— Initial Measurements in Fresh Water 

The ten-gallon aquarium was filled with fresh deionized 
water to conduct initial measurements and calibrations. Fresh water 
was used so that the possibility of electrical “cross-talk” between 
transducers would be reduced during these initial measurements. The 
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receive transducer was taken from its mount and used to test for 
electrical cross-talk by removing only the transducer face from the 
water. In this position, the receive transducer detected only a small 
portion of the imaging signal present in the water. By sliding the 
transducer in and out of the water, it was determined that the signal 
detected was not electrical cross-talk but was mostly noise picked up 
by the receiver through the sides of the transducer casing. 

The use of deionized fresh water did create a problem for 
bubble generation by electrolysis. The fresh water did not have 
enough conductivity to support electrolysis. Therefore, a small 
amount of table salt was added to increase the conductivity of the 
water. This proved adequate for the generation of small bubble clouds. 

Two different types of bubble clouds were investigated during 
this phase. The first bubble cloud was generated off a coil of brass 
wire which contained two different diameters of wire. Higher voltages 
were used and the bubble cloud was large and dense. Large bubble 
clouds appeared to make data collection difficult due to the large 
amount of bubble screening that occurred outside the sample volume. 
For this reason, a second, smaller bubble cloud was used. The second 
type of bubble cloud was a narrow stream of bubbles passing through 
the sample volume. It was created by a small, single-diameter piece of 
exposed brass wire at lower voltages. The bubble stream was still fairly 
dense, as shown in Figure 16, but did not appear to produce as much 
screening. These results will be discussed further in Chapter 5. 
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Figure 16. Narrow Bubble Stream in Aquarium 






The small, narrow bubble stream was used throughout the remainder 
of the experiment for all data collection. 

Transducer frequency response checks were done on both 
the imaging and pump transducers using both a calibrated Celesco 
LC-10 and a Naval Postgraduate School MA-1 hydrophone. The 
imaging transducer investigation was done to check for proper 
operation and to get an approximate idea of what frequency produced 
the highest imaging pressure amplitudes. The MA-1 was placed in 
the far field of the imaging transducer’s beam pattern, as shown in 
Figure 17. The frequency response was measured for imaging fre- 
quencies ranging from 2.20 MHz to 2.55 MHz at a far-field distance of 
6.7 centimeters. Imaging sound pressure amplitudes (pi) are calcu- 
lated in Appendix E from the measured sound pressure levels. The 
pump transducer frequency response was carefully checked at several 
distances using the LC-10 hydrophone as shown in Figure 18. The 
near field frequency response was recorded at one, two, four, and 
seven centimeters. The far field frequency response was recorded at 
ten and fourteen centimeters. The pump transducer frequency 
responses collected at various distances from the transducer face 
reflect the pump sound pressure amplitudes (pp) available. The pump 
sound pressure pp was calculated as shown in Appendix D. The best 
distance appeared to be a balance between the pump pressure 
amplitude available and the far field considerations discussed in 
Chapter II. Therefore, for most of the data runs, the sample volume 
was placed at a distance of ten centimeters from the pump transducer 
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Figure 17. Calibrated MA-1 Hydrophone For Measuring 
Imaging Transducer Frequency Response 
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face, which is in the far field for freqeuncies up to approximately 
300 kHz. 



A second check concerning the pump sound pressure versus 
frequency sweep time was also run. The LC-10 was placed ten 
centimeters from the -pump transducer face and the frequency 
response measured at' four different sweep times. A sizeable differ- 
ence in pump sound pressure levels was recorded between the 0.01- 
second sweep time and the higher sweep times. The most favorable 
sweep time appeared to be 0.10 seconds. The measurements for the 
various distances and sweep time frequency responses are also avail- 
able in Appendix D. Since the highest sound pressure levels (SPL) 
resulted with the sweep time of 0.10 seconds, that time was used 
throughout. 

The entire group of three transducers was placed in the 
aquarium such that standing waves and reflections produced a mini- 
mum amount of noise at the receive transducer. Absorptive material 
was placed in front of the transducers and at the surface of the water. 
This was done to help eliminate reflections and break up standing 
waves. An arrangement using absorptive material is shown in 
Figure 19. 

The placement of the pump, imaging, and receive trans- 
ducers relative to each other was investigated using three different 
arrangements. The first two arrangements were with the pump 
transducer in an upward-facing position (UP) and the sample volume 
at a distance of ten centimeters from the pump transducer face. The 
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Figure 19. Phase One Set-up With Absorptive Material 





sample volume was placed in either the near (Image Far Field or IFF) 
or far fields (Image Near Field or INF) of the imaging and receive 
transducers. The third transducer arrangement was with the pump 
transducer in a downward-facing position (DP). Here, the sample 
volume was put in only the far field of the imaging and receive 
transducers. Several sets of data were collected for each transducer 
arrangement. 

In all cases, the sample volume, bubble stream, and pump 
transducer sound fields for each arrangement were aligned using the 
LC-10 hydrophone. The pump transducer input level was 42 volts 
peak to peak and the imaging transducer input level was 6.0 volts 
(zero to peak), unless otherwise noted. 

Phase one data collection began with a fresh-water electrical 
cross-talk check. A pre-amplifier ground check provided a baseline 
noise level. The pump transducer was aimed directly into the face of 
the receive transducer to check for modulation and sidebands in the 
presence of no bubbles. A second “no-bubble test” to check for false 
sideband pressure amplitudes was conducted using an aluminum rod, 
or dowel, which acted as a reflector within the sample volume. The 
aluminum rod test is shown in Figure 20. Neither test produced dual- 
frequency sidebands. Finally, the receive transducer frequency 
response was recorded using bubbles as a reflector in the near field. 
This last check was used later to determine a more ideal image 
frequency which would avoid lower-frequency harmonics. 
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Figure 20. Aluminum Dowel Reflector Test Set-up 




The first fresh-water transducer arrangement used to record 
data was the Aquarium Upward Pump Image Far Fieldf, designated 
AQUPIFFf, where the subscript stands for fresh water. This arrange- 
ment placed the pump transducer face 10 centimeters from the sam- 
ple volume and maintained the far field for the high-frequency 
receiver. The following data runs were recorded using the AQUPIFFf: 

1. Wide Spectrum Anal 5 rzer Frequency Window— Pump and imaging 
sound fields for bubbles present versus no bubbles present. Pump 
sweep frequencies 15 kHz to 310 kHz. Imaging frequency 2.25 
MHz. 

2. Wide Spectrum Analyzer Frequency Window— Bubbles always 
present for imaging sound field only versus both pump and 
imaging sound fields. Pump sweep frequency 15 kHz to 310 kHz. 
Imaging frequency 2.25 MHz. 

3. Narrow Spectrum Anal 5 rzer Frequency Window— A “close-up” of 
the dual-frequency sidebands P+ and P_ with pump sweep fre- 
quencies of 15 kHz to 310 kHz. Imaging frequency 2.25 MHz. 

4. Repeat of run number three (3) with reduced pump power, i.e., 
less than 42 volts peak to peak for the pump power supply. 
Imaging frequency 2.25 MHz. Pump sweep frequencies of 15 
kHz to 310 kHz. 

5. Narrow Spectrum Analj^er Frequency Window— A “close-up” of 
the imaging sidebands P+ and P_ with an expanded pump sweep 
range of 15 kHz to 500 kHz. Imaging frequency 2.25 MHz. 

The second fresh-water transducer arrangement was the 
Aquarium Upward Pump Image Near Fieldf, designated AQUPINFf. 
This arrangement was used to record any increased signal sensitivity 
that could be gained by moving the receiver closer to the bubble cloud 
or stream. The pump transducer face was still maintained at a dis- 
tance of ten centimeters from the sample volume so the sample 
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volume would remain in the far field of the pump transducer beam 
pattern. The following data was recorded using the AQUPlNFft 

1. Narrow Spectrum Analyzer Frequency Window— Bubbles present 
and radiated by both pump and imaging sound fields. Pump 
sweeping frequencies of 15 kHz to 310 kHz. Imaging frequency 
2.25 MHz. 

2. Narrow Spectrum Anal3^er Frequency Window— Bubbles present 
and radiated by both pump and imaging sound fields. Reduced 
image level to 3.0 volts and expanded pump sweep frequencies of 
15 kHz to 420 kHz to study effects of harmonics. Imaging fre- 
quency 2.25 MHz. 

3. Repeat run two (2) except used receive transducer response data 
to move imaging frequency to more favorable 2.51 MHz to better 
avoid harmonics. Image level 3.0 volts. Pump sweep frequencies 
of 15 kHz to 420 kHz. 

The third and last fresh-water aquarium arrangement was the 
Aquarium Downward Pump Imaging Far Fieldf, designated AQDPIFFf. 

This arrangement was used to investigate any advantage gained by 
having the pump sound field oppose the bubble flow and by having the 
pump transducer face closer to the sample volume. A single set of 
data was recorded using the AQDPIFFf: 

Wide Spectrum Analyzer Frequency Window— Bubbles Present- 
Imaging sound field versus imaging plus pump sound fields, and the 
difference. Pump frequencies swept 15 kHz to 310 kHz. Imaging 
frequency 2.25 MHz. 

2. Phase Two— Seawater Testing 

The ten-gallon aquarium was filled with seawater. The 

electrical cross-talk checks conducted in phase one were repeated. 

The same procedure of immersing the transducer up to the 

transducer face and noting the response on the spectrum analyzer was 

used. The imaging transducer frequency response was again measured 
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in salt water using the MA-1 hydrophone. It did not differ from that 
measured in fresh water. 

Phase two data collection repeated that of phase one with the 
bubble sample volume in the imaging and receive transducer far field. 
The pump transducer was upward facing and the transducer face was 
at a distance of ten centimeters from the sample volume. In salt 
water, this configuration is designated AQUPIFFs- This data was 
collected for comparison with the fresh-water phase one data and, 
therefore, was run using the same variations. Data collected using the 
AQUPIFFs configuration includes: 

1. Wide Spectrum Analyzer Frequency Window— Bubbles present- 
imaging sound field only versus imaging plus pump sound fields, 
and the difference. Pump sweep frequencies of 15 kHz to 310 
kHz. Imaging frequency 2.25 MHz. 

2. Narrow Spectrum Analyzer Frequency Window— Single plot of 
bubbles radiated by both pump aind imaging sound fields. Pump 
sweep frequencies of 15 kHz to 310 kHz. Imaging frequency 
2.25 MHz. 

The aquarium filled with seawater was also used for the initial testing 
of the final transducer arrangement. This final arrangement was 
designed to be used in a large (3.0 meters deep) acoustically insulated 
tank. The Large Tank Downward Pump Image Far Fields, designated 
LTDPIFFs, is shown in Figure 15. The subscript “s” denotes the phase 
two use of salt water. The sample volume for this arrangement is 
approximately 7.0 centimeters from the pump transducer face. 
Appendix D also calculates the pump pressure amplitude (pp) for this 
distance. The pump pressure amplitude at 7.0 centimeters distance 
compares favorably with the pressure amplitude at 10 centimeters 
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distance. As mentioned previously, the imaging and receive 
transducers place the sample volume in the far field of the beam 
pattern. 



The following data runs were taken with the LTDPlFFs: 

1. Wide Spectrum Analyzer Frequency Window— Bubbles present- 
imaging sound field versus imaging plus pump sound fields, and 
the difference. Pump sweep frequencies of 15 kHz to 310 kHz. 
Imaging frequency 2.25 MHz. 

2. Wide Spectrum Analyzer Frequency Window— Bubbles present- 
imaging and pump sound fields. Imaging frequency shifted up to 
2.50 MHz to avoid harmonics. Pump sweep frequencies of 15 
kHz to 310 kHz. 

3. Wide Spectrum Analyzer Frequency Window— Bubbles present— 
Same set-up and format as run (1)— Extra imaging level to 9.83 
volts to bring out dual -frequency sidebands. Pump sweep 
frequencies of 15 kHz to 310 kHz. Imaging frequency 2.25 MHz. 

3. Phase Three— Large Acoustically Insulated Tank 

The third and final phase uses a large fresh-water tank, 
which is acoustically insulated, to investigate the possibility of reflec- 
tions causing excessive noise in the aquarium (phases one and two). 
The large tank provided an acoustically quiet environment due to its 
size, insulation, and the relatively low amplitude of the imaging and 
pump sound signals. The transducer arrangement was the same one 
used in the later portion of phase two. The Large Tank Downward 
Pump Image Far Fieldf (LTDPIFFfi, was used in the large tank with the 
same equipment rack as was used in the first two phases of the 
experiment. Figures 21 and 22 show the phase three experiment 
arrangements. 
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Figure 21. Large Tank Transducer Arrangement (LTDPIFFf) and Large Tank Access 




Figure 22. LTDPIFFf and Equipment Set-up Operational 
in Large Acoustically Insulated Tank 
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The LTDPIFFf data runs include the following: 

1. Wide Spectrum Analyzer Frequency Window— Bubbles present- 
imaging sound field only versus imaging plus pump sound fields, 
and the difference. Pump sweep frequencies of 15 kHz to 310 
kHz. Imaging frequency 2.25 MHz. 

2. Narrow Spectrum Analyzer Frequency Window— Sum Sideband 
(P+)— Bubbles present— Imaging sound field only versus imaging 
plus pump sound fields, and the difference Pump sweep fre- 
quencies of. 15 kHz to 310 kHz. 

The data and results for all data runs conducted in phases 
one, two, and three will be presented and briefly discussed in 
Chapter IV. 
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IV. RESULTS 



The results of the data collected during the course of the experi- 
ment are presented in this chapter. The figures that are the result of 
the data runs discussed in the Experimental Procedure section of 
Chapter III are plots of the HP3585A Spectrum Anal3^er CRT display. 
The actual figures were reproduced from the anal}^er by the HP 7090 
Measurement Plotting System. Each figure is a plot of sound pressure 
level (dB re 1.0 mW) versus frequency (MHz). The sound pressure 
levels can be easily converted to pressure amplitudes using the same 
methods described in Appendices D and E. The resolution bandwidth 
(RBW) from the CRT display is shown on each figure. 

This chapter is organized to present the results in the order that 
the data was taken during the experiment. A general discussion of 
each data run is included. The first data sets are for electrical cross- 
talk checks (both fresh and seawater), amplifier noise level measure- 
ments, measurement of the optimum receiver frequency response 
using bubbles as a reflector in the sample volume, a frequency 
modulation investigation using direct radiation of pump sound field 
energy, and comparison of the signals received by the receive trans- 
ducer for bubbles versus a solid reflecting object in the sample volume. 
The data runs from phases one, two, and three of the experimental 
procedure follow the discussion of the initial calibrations. 

Electrical cross-talk checks were conducted in both fresh water 
and seawater to determine if the increased conductivity of seawater 
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would cause a portion of the electrical Signal to be transmitted from 
the imaging transducer to the receive transducer. Figures 23a and 
23b are measurements of the signals received via the casing of the 
receive transducer at 2.25 MHz. Both figures are in terms of sound 
pressure levels (SPLs). Both figures are the difference between the 
2.25 MHz signal being present and the noise when no signal is pres- 
ent. The signal spike at 2.25 MHz is approximately 16 dBm for both 
fresh water and seawater. It is concluded, therefore, that electric<»l 
cross-talk must be minimal and is the same for both fresh water and 
seawater. A more likely source of the signal spike is energy that is 
transmitted acoustically through the receive transducer casing. 

A baseline noise measurement was performed by grounding the pre- 
amplifier input. The results are plotted in Figure 24. The noise level from 
150 kHz up to 3.50 MHz is constant at approximately -103 dB re 1.0 mW. 

The receive transducer frequency response using swept frequen- 
cies of the sound field reflected from a bubble stream in the sample 
volume is plotted in Figure 25. Figure 25 was used to determine the 
adjustment of the imaging frequency (fi) from 2.25 MHz to 2.50 MHz 
in order to better avoid harmonics of the pump sound field and 
resonating bubbles. Figure 25 shows that the receive transducer 
frequency response begins to fall off after 2.50 MHz. 

A frequency modulation check in fresh water was conducted to 
investigate the effects of direct radiation of the pump sound field on 
the face of the receive transducer. The pump transducer was pointed 
directly at the receive transducer and the difference between the 
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Figure 23a. Fresh Water Electrical Cross-Talk Check at Imaging Frequency 
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Figure 24. Baseline Noise Measurement 
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' Frequency (MHz) RBW: 10.0 

Figure 25. Receive Transducer Frequency Response Due to 

Reflected Imaging Sound Fields (Swept Frequencies) 



signal and noise plotted in Figure 26. This figure demonstrates that har- 
monics of the pump (fp = 15 kHz to 310 kHz) exist throughout the fre- 
quency spectrum and higher than the imaging frequency of 2.25 MHz. 
The pump transducer was driven at a maximum input level of 42 volts 
peak to peak. This did produce very low pressure levels at the side- 
bands about the imaging frequency, but they were considered insig- 
nificant to bubble detection because it required direct pump sound field 
radiation to cause modulation and produce those low-level sidebands. 

The final calibration check was to measure the sound pressure 
levels of the imaging and pump sound fields as they were reflected 
from a solid aluminum dowel placed in the sample volume. The 
purpose was to see if the combined sound fields would modulate if 
reflected from a solid. Figure 27a is the noise measurement of only 
the imaging sound field as it is reflected to the receive transducer. 
The amount of sound energy reflected at the peak (2.25 MHz) was 
adjusted to be similar to the same level as that reflected by a bubble 
stream (SPL approximately -27 dB re 1.0 mW). Figure 27b is a plot of 
both the imaging and pump sound fields reflected by the dowel. 
Figure 27c is a plot of the difference between the reflected dual sound 
field and the imaging “noise" sound field. Figure 27c shows no 
sidebands or harmonics created by the aluminum dowel, which is as 
expected since the dowel is solid and does not resonate. 

All of the previously mentioned initial measurements were con- 
ducted to prove that the signals received in the data sets of phases 
one, two, and three were caused by bubbles that produced legitimate 
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Figure 26. Frequency Modulation Due to Direct Radiation of Pump Sound Field 
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Figure 27a. Sideband and Harmonic Investigation Using Solid Reflector— Imaging Sound Fie 
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sidebands about the imaging frequencies. The following sections pre- 
sent the data as it was taken during the experimental procedure and 
provide some observations. The same abbreviations for transducer 
arrangements used in Chapter III apply here in Chapter IV. 

A. PHASE ONE DATA 

Figures 28a, 28b, and 28c represent a baseline measurement 
which compares the received imaging, pump, and resonating bubble 
sound signals in the presence of bubbles versus the absence of bubbles. 
These figures show data for the Aquarium Upward Looking Pump 
Image Far Fieddf (AQUPIFFf), where the subscript “F is for fresh 
water, set-up and a wide HP3585A spectrum analyzer frequency 
window (0.0 MHz to 2.75 MHz). Figure 28a shows sound pressure 
levels for the pump sound field being swept from fp = 15 kHz to 310 
kHz and a sharp imaging sound level “spike” at 2.25 MHz. No bubbles 
are present in the sample volume and, therefore, these signals are 
considered “noise." Figure 28b demonstrates the change when 
bubbles are introduced into the sample volume. Harmonics are now 
evident in the frequency spectrum almost up to the imaging frequency 
(fi). Sharp sideband levels appear on either side of the imaging 
frequency (fi) and extend up and down the spectrum to fi ±310 kHz. 
It can be noted here that the signal-to-noise ratio for the imaging 
sound pressure level is approximately 70 dB re 1.0 mW above the 
“noise" level. Figure 28c is the signal difference between the 
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Figure 28a. Wide Spectrum Analyzer Window— AQUPIFFf— 

Imaging and Pump Sound Fields— No Bubbles 
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Figure 28b. Wide Spectrum Analyzer Window— 'AQUPIFI^— 

Imaging and Pump Sound— Bubbles Present 
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Figure 28c. Wide Spectrum Analyzer Window— AQUPIFFf— 

Bubble/No Bubbles Difference 





the presence of bubbles (Figure 28b) and the absence of bubbles 
(Figure 28a). 

Figure set 29 is the result of constantly maintaining the bubbles as 
reflectors in the sample volume and shows the difference between the 
presence and absence of the pump sound field while the bubbles are 
irradiated by the imaging sound field. Data was collected in this 
manner for the duration of the experiments. The objective was to 
provide a constant target in the sample volume and maintain a 
constant peak imaging sound pressure level signal -to-noise ratio. 
Figure 29a is a “noise” measurement of the imaging sound field as it is 
reflected off of the bubbles in the sample volume. In Figure 29b, the 
pump sound field had been turned on and created distinct harmonics 

and dual-frequency sideband pressure levels (that can be converted to 
P+ and PJ. Note again that the harmonics go high up the frequency 

spectrum to the point of the imaging frequency (fj, which suggests 
that increasing the imaging frequency above 2.25 MHz would help 
avoid the harmonics. Figure 29c is a plot of the difference in signal for 
the presence of both imaging and pump sound fields versus only the 
imaging sound field. 

Figure 30 is a narrow spectrum analyzer window with both 

imaging and pump sound fields present, creating the dual-frequency 
sideband sound pressure levels (P+ and PJ about the imaging fre- 
quency pressure level. Figure 30 shows a sizeable signal-to-noise ratio 
of approximately 12 to 15 dB for these sidebands (at ± fp). The 
sideband pressure levels have sharp cut-offs on either side of the 
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Figure 29a. Wide Spectrum Analyzer Window— AQUPIFFf— Imaging Sound Field Only 
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Figure 29b. Wide Spectrum Analyzer Window— AQUPIF^ —Imaging and Pump Sound Fields 




Ov^ui O’ I 3-1 SP) I3A3q 9jnss9Jj punos 

115 



Figure 29c. Wide Spectrum Analyzer Window— AQUPIFFf 
Imaging and Pump/Imaging Only Difference 
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Figure 30. Narrow Spectrum Analyzer Window — AQUPIFIJ — Imaging and Pump Sound Fields 





imaging pressure level “spike” at fi ± fp where the pump sound field 
comes to the end of its sweep (fp = 310 kHz in this case). This sharp 
cutoff indicates the presence of even smaller bubbles than those 
resonated at 310 kHz. 

Figure 31 represents the same data as Figure 30, except that the 
pump frequency sweep range was expanded so that the pump sound 
field was swept from fp = 15 kHz to 500 kHz at a normal level of 42 
volts peak to peak. This was done to study the effects of pump sound 
field harmonics. The previous figures show sharp sideband cutoffs at fi 
±310 kHz, indicating that bubble sizes exist beyond the 310 kHz reso- 
nance range. Figure 31 shows the effect, or trade-off, of the expanded 
pump sweep frequencies (fp) since the harmonics carry beyond the 
imaging frequency (fi) and cover the sideband (fi ± fp) pressure levels 
with noise. 

The possibility of reducing the pump sound field harmonics dis- 
cussed in the previous paragraph exists at reduced pump sound field 
power. Figures 32a and 32b are the plots of the results when reducing 
the pump transducer input from 42 volts peak to peak to 20 volts 
(Reduced Pump Power (1)) and 10 volts (Reduced Pump Power (2)), 
respectively. This power reduction lowered the pump pressure 
amplitude (pp) supplied to the bubbles and, therefore, reduced the 
harmonics. Reducing the pump pressure amplitudes also lowered the 
pressure amplitudes of the dual-frequency sideband. The expanded 
pump sweep frequencies of fp = 15 kHz to 500 kHz were maintained 
for the plotting of both Figures 32a and 32b. 
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Figure 31. Narrow Spectnmi Analyzer Window— AQUPIFFf 
Imaging and Expanded Pump Soimd Field 
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Figure 32a. Narrow Spectrum Analyzer Window — AQUPIFI^ — Imaging and 
Eicpanded Pump Sound Fields — Reduced Pump Power (1) 
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Figure 32b. Narrow Spectrum Analyzer Window — AQUPIFIJ — Imaging and 
Expanded Pump Sound Fields— Reduced Pump Power (2) 



The transducer arrangement was changed to the Aquarium 
Upward Looking Pump Image Near Field in fresh water, designated 
the AQUPINFf. This was done to investigate the receive transducer 
sensitivity and the imaging transducer power inputs when the sample 
volymie was in the near field. Figure 33 is a narrow spectrum analyzer 
window plot of the sideband pressure levels and associated noise 
levels for the normal imaging level of 6 volts and imaging frequency 
(fi) of 2.25 MHz. The pump frequency (fp) was swept from 15 kHz to 
310 kHz at 42 volts peak to peak. The sideband pressxire levels shown 
in Figure 33 show a clear cut-off at fi ± 310 kHz, indicating the 
termination of the pump frequency sweep, not the lack of smaller 
bubbles in the sample volume. 

Since the bubble sizes present in the sample volume extend 
beyond the 310 kHz resonance frequency range. Figure 34 is plotted 
for the AQUPINFf with the input to the imaging transducer reduced to 
3.0 volts and the pump frequencies (fp) expanded to range from 15 
kHz to 420 kHz. The imaging input level was reduced to check the 
effects on the harmonics created by the resonating bubbles and 
demonstrate the direct relationship of the imaging sound field (pi) to 
the sidebands (P+). Figure 34 shows that the sidebands (fi ± fp) are 
still present at the reduced image transducer input level, however, the 
resonance harmonics still go higher up the frequency spectrum 
beyond the upper sideband (fi ± fp) pressure level at 2.67 MHz. Note 
that reducing the imaging pressure amplitudes lowers the dual- 
frequency sideband pressure amplitudes (see Figure 33 for comparison). 
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Figure 33. Narrow Spectrum Analyzer Window — AQUPINFj — Imaging and Pump Sound Fiel 
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Figure 34. Narrow Spectrum Analyzer Window— AQUPINFf— Imaging and 
Eicpanded Pump Soimd Field — Reduced Image Power 





In an effort to avoid the noise caused by harmonics, the imaging 
frequency (fi) was increased to 2.51 MHz, moving the imaging 
pressure level peak further up the frequency spectrum. The results of 
this frequency increase are plotted in Figure 35. The imaging input 
level is still reduced to 3.0 volts and the pump frequencies (fp) are still 
expanded to range from 15 kHz to 420 kHz so that a comparison with 
Figure 34 can be made. Note in Figure 35 that the upper sideband 
pressure levels (at fi + fp = 2.67 MHz) are now free of the resonating 
bubble harmonic noise. 

The final data set collected in phase one was for the Aquarium 
Downward Looking Pump Image F^ Field in fresh water, designated 
AQDPIFFf. This data set is plotted in Figures 36a, 36b, and 36c. The 
downward-looking pump transducer experiment was for comparison 
with the upward -looking pump transducer experiment. The trans- 
ducer inputs are: imaging transducer at fi = 2.50 MHz and input level 
at 6.0 volts; pump transducer at fp = 15 kHz to 310 kHz and input 
level at 42 volts peak to peak. The imaging frequency (fi) is high on 
the frequency spectrum to better avoid the harmonics induced by 
resonating bubbles and the pump sound field. 

Figure 36a is a wide spectrum analyzer window plot of the imag- 
ing sound field reflected from a small stream of bubbles in the sample 
volume. Figure 36a is the “noise” level for the downward-looking 
pump data runs. Figure 36b is a plot of the imaging and pump sound 
fields interacting with the bubble stream. Figure 36c is the difference 
between the signal levels plotted in Figures 36b and 36a. It appears 
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Figure 36a. Wide Spectrum Analyzer Window— AQDPIFIJ— High Imaging Sound Field Only 
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Figure 36b. Wide Spectrum Analyzer Window — AQDPIF]^ — High Imaging and Pump Sound Fields 
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from the. pressure levels plotted in Figures 36b and 36c that the 
downward-looking pump transducer provides the same signal-to- 
noise ratios as the upward-looking pump transducer. 

B. PHASE TWO DATA 

Phase two was conducted in the ten-gallon aquarium with sea- 
water and two different transducer arrangements. Further transducer 
arrangements were not used in this phase because the primary pur- 
pose of phase two was only to show the compatibility of the dual- 
frequency bubble density measurement technique with seawater and 
the ship wake bubble density problem. To draw a comparison between 
fresh water and seawater, the Aquarium Upward Looking Pump Image 
Far Fields (AQUPlFFs) (the subscript “s” designates each experiment 
conducted in seawater) transducer arrangement was used to repeat 
the previous fresh-water data runs. The second transducer 
arrangement used in seawater was the Large Tank Downward Looking 
Pump Image Far Fields (LTDPIFFs) mount designed for use in the 
large, acoustically insulated tank (approximately three meters deep). 
This second transducer arrangement was initially tested in seawater 
and, because it closely approximated the fresh water downward- 
looking pump transducer mount (AQDPIFFd, was used to compare the 
downward-looking pump in seawater versus fresh water. 

The first seawater data plotted was obtained using the AQUPlFFs 
transducer arrangement. A wide spectrum analyzer window, with 
standard imaging and pump transducer inputs (fi = 2.25 MHz at 6.0 
volts, fp = 15 kHz to 310 kHz at 42 volts peak to peak), was used to 
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plot Figures 37a, 37b, and 37c. As before. Figure 37a was plotted with 
only the imaging sound field reflected from the bubble stream, and is 
considered the “noise” level. Figure 37b plots the received signal result- 
ing from bubbles being irradiated by both imaging and pump sound 
fields. Figure 37c is the difference between the signal received from the 
bubbles with both sound fields present (Figure 37b) ^d the imaging 
noise (Figure 37a). The sideband pressure levels present at the side- 
band frequencies (fi ± fp) correspond to the pump frequencies that 
were swept, however, there appears to be a lower signal-to-noise ratio 
for both imaging and sideband pressure levels in seawater. Figure 38 
is a plot using a narrow spectrum analyzer window, or “close-up,” of 
the sideband pressure levels shown in Figure 37b. Figure 38 data also 
shows the slightly lower imaging and sideband signal-to-noise ratio 
when compared to the same data of Figure 30 for fresh water. 

The second transducer arrangement used in phase two of the 
experiment was the LTDPIFFs in the aquarium. Rgures 39a, 39b, and 
39c again show harmonics which result from resonating bubbles under 
the influence of the pump -sound field (sweeping frequencies fp = 1 5 
kHz to 310 kHz), and the sideband pressure levels which result from 
dual-frequency bubble excitation. Figure 39a is a plot of the “noise” 
signal received when the bubbles are irradiated by only the Imaging 
sound field. Figure 39b is a plot of the received signal from the 
bubbles when both imaging and pump sound fields are present. Figure 
39c is a plot of the reflected imaging “noise” subtracted from the 
imaging and resonating bubble sound fields. Figure 40 is a 



130 




0*1 9P) l3A9q 3jnss3Jd punos 



131 



Figure 37a. Wide Spectrum Analyzer Window— AQXJPIFT's~ Imaging Soimd Field Only 



CO 




lO 

I-H 

I 



lO 



(^ui 0*1 aP) pA9q 9jnss9Jd punos 

132 



I 



Figure 37b. Wide Spectrum Analyzer Wndow—AQUPIFI^— Imaging and Pump Sound Fields 
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Figure 38. Narrow Spectrum Analyzer IX^dow—AQUPIF^ — Imaging and Pump Sound Fields 
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Figure 39a. Wide Spectrum Analyzer Window— LTDPIFFs— Imaging Sound Field Only 
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Figure 39b. Wide Spectrum Analyzer Window— LTDPIFFs— Imaging and Pump Sound Fields 
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Figure 39c. Wide Spectrum Analyzer Window— LTDPIFFs 
Imaging and Piunp/Imaging Only Difference 




<N 



lO 

I 



in 

7 



lO 

t 



lO 

1 



lO 



N 

I 



00 

I 



lO 

o 

I 



o 



ic 

cs 

I 



0‘t 9P) pA^i 9jnss9Jd punos 

138 



I 



Figure 40. Wide Spectrum Analyzer \»^dow— LTDPIFFs— High Imaging and Pump Sound Fields 



wide -spectrum analyzer window repeat of Figure 39b, except the plot 
shows a change in the imaging frequency to 2.50 MHz so as to move 
the dual-frequency sidebands away from the resonance harmonics. 
The plots in Figures 39a, 39b, and 39c were all recorded using an 
imaging frequency at 2.25 MHz at 6.0 volts of transducer input, and 
swept pump frequencies of 15 kHz to 310 kHz set at 42 volts peak to 
peak transducer input. Again, the dual-frequency sideband sound 
pressure levels in both Figures 39b and 40, while both at 
approximately -90 dB re 1.0 mW, appear to be slightly lower than the 
leyels recorded in fresh water. The sideband pressure levels plotted 
in Figure 40 do have a better signal-to-noise ratio than those plotted 
in Figure 39b due to the higher imaging frequency used. The higher 
imaging frequency avoids the resonance harmonics of the bubbles. 

C. PHASE THREE DATA 

Phase three was conducted in one of the Naval Postgraduate 
School’s large, fresh water, acoustically insulated tanks. The object of 
this phase was to Investigate the effects of a reflection-free environ- 
ment on the dual-frequency sideband sound pressure level signal-to- 
noise ratios. The Large Tank Downward Pump Image Far Fieldf in 
fresh water was used in the non-reflecting tank. 

The first series of data is plotted in Figures 41a, 41b, and 41c. 
The same data format of phases one and two is repeated here. Figure 
41a is a plot of the received reflected imaging signal from the bubbles 
and is considered “noise.” Figure 41b is the received signal from the 
bubbles with both imaging and pump sound fields present. Figure 41c 
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Figure 41a. Wide Spectrum Analyzer Window— LTDPIFFf— Imaging Sound Field Only 
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Figure 41b. Wide Spectrum Analyzer Window— LTDPIF^ —Imaging and Pump Sound Fields 
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Figure 41c. Wide Spectrum Analyzer Window— LTDPIFFj — 

Imaging and Pump/Imaging Only Difference 




is the difference between the signal from the bubbles in the presence 
of imaging and pump soimd fields and the imaging “noise” only. 

■ The second series of phase three data, plotted in Figures 42a, 
42b, and 42c, is a narrow spectrum analyzer window, or “close-up,” of 
the first series. Figure 42a is a plot of the imaging noise. Figure 42b 
is plotted with both imaging and pump sound fields radiating the bub- 
bles, and Figure 42c is the difference between the latter two. All 
phase three experiments utilized an imaging frequency of 2.25 MHz at 
a transducer input level of 6:0 volts. The pump frequencies were 
swept from 15 kHz to 310 kHz. The pump transducer input level was 
42 volts peak to peak. 

The results from all three experimental phases has been briefly 
discussed in this chapter. Chapter V, Conclusions and Recommenda- 
tions, will use these results and further discuss the implications of the 
data in a more general sense. 
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2.250 2.300 2.350 2.400 2.450 2.500 2.550 2.600 2.650 2.700 2.750 

Frequency (MHz) RBW: 1.0 KHz 

Figure 42a. Narrow Spectrum Analyzer Window— LTDPIFFj— 

Upper Sideband— Ima^g Sound Field Only 
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2.250 2.300 2.350 2.400 2.450 2.500 2.550 2.600 2.650 2.700 i 

Frequency (MHz) RBW: 1.0 

Figure 42b. Narrow Spectrum Analyzer Window — LTDPIFFj — 

Upper Sideband— Imaging and Pump Sound Fields 
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2.250 2.300 2.350 2.400 2.450 2.500 2.550 2.600 2.650 2.700 2.750 

Frequency (MHz) RBW: 3.0 KHz 

Figure 42c. Narrow Spectrum Analyzer Window— LTDPIFIJ— Upper Sideband- 
Imaging and Pump/Imaging Only Difference 



V. CONCLUSIONS AND RECOMMENDATIONS 



The goal of this thesis is to provide a theoretical and experimental 
investigation into the feasibility of using the Dual Frequency Pump 
Technique to acoustically determine point-by-point bubble density 
distributions within surface ship wakes. The theoretical investigation 
reviewed bubble resonance, dual-frequency bubble excitation, and 
transducer beam patterns in Chapter II. The experimental considera- 
tions, such as transducer orientation, sample volume placement, bub- 
ble screening, and idealized transducer inputs were discussed in 
Chapter III. Data that was gathered during the experiment and some 
general observations concerning that data are presented in Chapter IV. 
Some general conclusions and recommendations are presented here 
in Chapter V. 

A. CONCLUSIONS 

The Dual Frequency Pump Method of acoustically detecting bub- 
bles and determining point-by-point bubble cloud densities is very 
promising for surface ship wake measurements. The design of an 
actual measurement device requires a far more intensive and detailed 
research effort, but the dual frequency technique is practical for sur- 
face ship wake bubble density measurements. This thesis demon- 
strates that practicality through the discussion of the following topics. 

1. Bubble Resonance 

Bubbles undergoing resonance and volume pulsations are the 
keys to making the dual frequency pump technique work. Bubbles 
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resonate easily in both fresh and salt water as is evidenced by the 
presence of the harmonics of the resonance sound field plotted in 
Figures 28b and 37h. The harmonics verify that the bubbles are acting 
as a sound source. The data collected during the seawater experi- 
. ments in the aquarium shows that dual frequency excitation works in 
seawater; but the data plotted in Figures 37 through 40 for seawater 
shows a slightly lower signal-to-noise ratios as discussed in the 
Chapter IV section on phase two data. The reason for a lower signal 
level in seawater is not clear. 

Figures 27a, 27b, and 27c demonstrate the Dual Frequency 
Pump Method’s ability to detect only bubbles, not solids, at the imag- 
ing frequency sidebands (coi ± cop). These figures are plots of the 
received spectrum when an aluminum dowel was used as a reflector 
instead of bubbles. Figure 27c is a plot of the sound pressure level of 
the reflected sound field for pump frequencies swept from 15 to 310 
kHz. The lower frequency spectrum of the plot shows the reflected 
sound signal but the upper spectrum shows no sideband signals, indi- 
cating correctly that there are no bubbles present. Seawater, particu- 
larly in ship wakes, is rich in suspended solids and biological life. The 
Dual Frequency Pump Method can distinguish between a bubble and a 
suspended solid. 

The solid reflector test was carried a step ftirther by radiating 
the face of the receive transducer directly with the pump sound field. 
Figure 26 plots no significant frequency modulation, which could give 
a false bubble indication, caused by the maximum pump sound 
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pressure amplitude in the absence of bubbles. The receive transducer 
was radiated with the pump sound field in the presence of an imaging 
signal. However, it should be noted that the imaging sound field, 
which acts as the “carrier” of the pump sound field, was not aimed 
directly at the receiver. These two simple tests help demonstrate that 
the received signals plotted in all the figures are in fact a product of 
bubble resonance and dual-frequency excitation. 

The results of Chapter 11 equations (7) and (8), which are 
presented in Tables la, lb, and Ic, deserve comment. First, in direct 
regard to dual-frequency bubble excitation,. these tables show the need 
for sweeping a wide range of pump frequencies. Shankar and 
Newhouse used only a limited pump sound field frequency range for 
their biomedical-based research [Refs. 3 and 12]. The arbitrary range 
of bubble radii selected for our demonstration of the surface ship wake 
measurement problem suggests a broad sweep range of approximately 
20 kHz to 500 kHz in order io compensate for increased hydrostatic 
pressures. Second, the hydrostatic pressure term in equations (7) and 
(8) is significant to surface ship wakes in that it could help explain 
why surface ship wakes persist longer than expected at lower fre- 
quencies. Tables la, lb, and Ic demonstrate that a bubble of given size 
has a much lower resonance frequency^ at the ocean surface where the 
hydrostatic pressure is low than at some depth where the pressure is 
greater. Remembering that the ship wake is defined by the presence 
of many different-sized bubbles of different resonance frequencies, the 
acoustic presence of the wake is determined by the large amounts of 
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absorption or scattering that occurs' at those frequencies. Logic would 
dictate that the lower-frequency bubbles, which are larger in radius, 
would disappear faster, either by rising to the surface more quickly or 
by breaking up into smaller bubbles of higher resonance frequency. 
The small, high-resonance frequency bubbles with the slow rise times 
should be present for longer periods. Wake profile data, which is 
taken in vertical “slices" or columns, supports the logic by showing 
the presence of high-frequency wakes. However, the wake data also 
shows the presence of a lower-frequency wake which persists much 
longer than expected. A possible explanation could be that, if a bubble 
of particular radius started its rise toward the surface from deep in the 
wake, the bubble’s resonance frequency becomes less as the bubble 
rises. A vertically oriented wake profiling device would not sense a 
change in the bubble’s position within the water column, only the 
presence of the lower resonance frequencies. The bubble rise mecha- 
nism is coupled to njany other variables such as turbulence, diffusion, 
and expansion of the bubbles, but it may offer a simplistic insight into 
a possible cause of low-frequency wake persistence. The Dual Fre- 
quency Pump Method of acoustic wake measurement could solve the 
lack of vertical position data by taking point measurements within the 
wake. 

Bubble resonance can also have an adverse effect on the Dual 
frequency Pump Method if the resonance creates strong harmonics. 
Excessive pump pressure amplitudes at the sample volume location 
results in resonance harmonics far up the frequency spectrum. If the 
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harmonics extend to high enough frequencies, they raise the noise 
level about the imaging frequency, which creates difficulties in 
detecting the low-amplitude sideband pressures. These resonance 
harmonics are recorded in most of the plots and discussed again in a 
later section. 

2. Bubble Screening and Sample Volume Size 

Early in the discussion of the experimental procedure, men- 
tion was made of the use of narrow bubble streams passing through the 
sample volume versus the generation of large bubble clouds. This was 
necessary due to the large amounts of bubble screening which take 
place in the presence of a large cloud. 

The sample volume created by the intersection of the imaging 
and receive transducer beam patterns is small relative to a large bub- 
ble cloud. A large bubble cloud will contain the sample volume at some 
point on the cloud’s interior. If this is the case, many bubbles are out- 
side the sample volume. This has two degrading effects. First, the 
imaging and pump sound undergo large amounts of attenuation due to 
screening bubbles prior to radiating the target bubbles in the sample 
volume. Second, with the resonating target bubbles acting as a sound 
source, the receive transducer has difficulty detecting the radiated 
bubble sound field due to attenuation from screening bubbles. The 
effect of screening bubbles is, therefore, twofold when the sample vol- 
ume is inside a large bubble cloud. This large cloud “mechanism” was 
supported by the results of initial data collection efforts. The receive 
transducer and spectrum analyzer were unable to detect any dual- 
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frequency sideband pressure levels while only barely detecting 
primary pump and imaging sound fields. 

The use of narrow bubble streams passing thurough the sample 
volume 5 delded much greater success. Dual-frequency sideband pres- 
sure levels were detected during all three phases of the experiment 
when the small stream was used. The narrow bubble stream success- 
fully eliminated the screening bubbles outside the sample volume. 
This provided all the sound fields clear paths to and from the trans- 
ducer faces, thus keeping attenuation to a minimum. This evidence 
emphasizes the need to keep the bubble streams to be measured small 
and contained within the sample volume. If large bubble clouds, such 
as ship wakes, are to be measured, it will be necessary to isolate the 
sample volume from the bubble cloud so that sound paths will not be 
interrupted by screening bubbles.. 

3. Dual Frequency Sideband Pressure Amplitudes 

Data collected in phases one, two, and three of our experi- 
ments confirms the finding of Shankar and Newhouse concerning the 
dual-frequency sideband pressure levels. Shankar and Newhouse 
found that, in order to detect bubbles at the dual-frequency sidebands 
(coi ± cop), it was necessary to maintain the ratio of the imaging signal 
level (equivalent to Pi) to sideband level (equivalent to P+) ratio of 
approximately 63 dB. Figures 29a, 29b, and 29c provide a graphic 
illustration of this signal level ratio. The imaging (Pi) signal level in 

Figure 29b is approximately 80 dB above the background noise and the 
dual-frequency sideband levels (P+ and PJ are approximately 20 dB 
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above ' the noise at their peak level. The • signal level difference 
between the two is approximately 60 dB. This ratio became important 
because, for every experiment that was conducted, the imaging sound 
level (Pi) needed to be at least 55 to 60 dB above the noise level 
before dual -frequency sideband levels could be detected. The side- 
band pressure amplitudes are very small compared to those of the 
reflected imaging and pump sound fields, therefore, a large signal-to- 
noise ratio is needed to make this technique successful. 

Table III theoretical results of equations (28a) and (28b) also 
support the conclusion that the signal-to-noise ratios of the imaging, 
dual-frequency sideband, and background pressure amplitudes are key 
factors in the success of the Dual Frequency Pump Method. Table III 
results are calculated assuming that the pump and image sound fields 
are of constant pressure amplitudes for all frequencies. This means 
that Pi and pp are constants in equations (28a) and (28b). The com- 
puted sideband pressure amplitudes compare favorably with those 
found by Shankar and Newhouse. Table III shows the sideband pres- 
sure amplitudes to be very small compared to the pump pressure 
amplitudes. The table also highlights the logical conclusions that both 
the sideband pressure amplitudes for a bubble at resonance and the 
noise pressure amplitude of the bubble not at resonance increase with 
bubble size. 

The data plotted for all figures was collected using the “max 
hold” feature of the HP3585A spectrum analyzer. The “max hold” 
feature stores the peak pressure level at each frequency on the CRT of 
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the analyzer. Therefore, all the plots represent bubble size distribu- 
tions rather than actual bubble counts. The pressure level peaks at the 
sidebands indicate that at least one bubble of that size was present 
during the 132 -second statistical sample time. 

The data plotted for all these phases of the dual-frequency 
experiment shows an increase in pressure level as the sideband fre- 
quencies increase and decrease about the center (carrier) imaging 
frequency. Figures 42a, 42b, and 42c (plots of the upper sideband) 
are good examples of the maximum sideband level being located at a 
corresponding resonance frequency of approximately 300 kHz (coi + 
cop = 2.55 MHz). This indicates the presence of a small bubble and 
appears contrary to the results of equations (28), in that small bubbles 
radiate small pressure amplitudes. The experimental results are not 
contrary, however, due to the fact that the actual pump pressure 
amplitude increases with frequency because of the pump transducer 
response. The experiment provided the smaller bubble with a greater 
pp than the theoretical calculation. The theoretical calculation was 
performed utilizing the pump pressure amplitude as a constant. Theo- 
retical results similar to those plotted in the figures could be obtained 
if the input (pp) to equations (28a) and (28b) is varied with frequency. 

The difficulties associated with obtaining actual bubble counts 
and maintaining a constant acoustic pump pressure level (to make 
theoretical calculation easier) can both be solved using software. A 
computer could individually store the dual-frequency sideband pres- 
sure levels on every sweep of the analyzer and add them over the 
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course of the statistical sample time. The computer could also com- 
pute theoretical pressure levels (P+) for a variety of pressure level 
inputs (pi and pp) and compare those results to the signals received 
from the sample volume at discrete frequencies. 

4. Imaging and Pump Sound Fields 

The importance of the pump sound field frequencies has 
already been discussed, however, the experimental results further 
show that the pump sound field amplitude, the imaging frequency, and 
the imaging sound field amplitudes are important. The imaging sound 
field amplitude has also already been mentioned in the previous sec- 
tion. The effect of two different imaging pressure levels can be 
observed in Figures 33 and 34, where the imaging input level was 
reduced for data in the Figure 34 plot. Reducing the imaging input 
voltage by half reduced the sideband levels by approximately 6 dBm 
and produced a less -favorable signal-to-noise level difference. This is 
as expected since equations (28) show that the upper sideband level 
(P+) is directly related to both imaging and pump inputs. 

The effects of the pump pressure amplitude on sideband 
pressure amplitude are plotted in Figures 31, 32a, and 32b. The 
phase one AQUPIFFf experimental data was recorded to demonstrate 
the effect reduced pump pressure amplitude had on harmonic and 
sideband pressure levels. High pump pressure amplitudes at the 42- 
volt peak-to-peak transducer input levels create harmonics which 
exist at high frequencies on the frequency spectrum. These harmon- 
ics create a noise level problem at the dual-frequency sideband (coi ± 



155 



(Op) locations. The harmonics, caused by an expanded pump frequency 
range (cop = 15 to 500 kHz) are also detected at both upper and lower 
dual-frequency sideband locations on the frequency spectrum (as 
plotted in Figure 31). The sideband pressure levels are barely above 
the noise of the harmonics. Therefore, reduced pump pressure 
amplitudes were tested and the results shown in Figures 32a and 32b. 
Figures 32a and 32b show a desirable reduction in harmonic noise. 
The dual-frequency sideband levels drop approximately 6 dBm, as 
expected (equation (28)), with each input level reduction. 

The ideal imaging frequency was also investigated in both 
fresh water and seawater. Several transducer arrangements were used 
to compare results. The data of Figures 35, 36a, 36b, 36c, and 40 
indicates that an imaging frequency higher than 2.25 MHz avoided the 
noise level that was created by harmonics. Increasing the imaging 
frequency was effective in avoiding harmonic noise from all transducer 
arrangements and both fluids. Surface ship wakes should have a low 
noise level at ultrasonic frequencies. Harmonics aire one of few effects 
that can raise the ultrasonic noise level. 

5. Sample Volume Location. Transducers, and Transducer 

Arrangements 

The sample volume should be located in the far field of all 
transducer beam patterns. The increase in sound field pressure 
amplitudes is not significant to justify moving the sample volume to 
the near field. These plots of data collected with the sample volume 
in the near field (AQUPINFf) of the imaging and receive transducers 
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(see Figures 33, 34, and 35) did not show an improvement over the 
experiments with the sample volume in the far field. The pump 
transducer was placed whenever possible so that the sample volume 
would be in the far field at the pump transducer. The LTDPlFFf and 
the AQDPlFFf arrangements required the pump transducer to be 
closer than the calculated far field point to the sample volume due to 
physical transducer arrangement limitations. Near field points were 
chosen carefully to provide equitable sound pressure levels to those of 
the far field. 

The transducers used for these experiments were low cost 
and did not have ideal center frequencies or frequency responses. 
This was particularly true of the pump transducer. A lower center 
frequency would have allowed a better frequency response about reso- 
nance frequencies of interest. The receive transducer (center fre- 
quency 2.25 MHz) was not -calibrated and therefore the pressure 
amplitudes of the dual-frequency sidebands could not be calculated 
from the HP3585A sound pressure level data. Calibrated receive 
transducers must be used to obtain sound pressure amplitudes. 

Several transducer arrangements were studied in phases one, 
two, and three. The largest physical difference between any one 
arrangement was the upward-looking and downward-looking pump 
transducers. No appreciable system performance difference was 
detected between the upward and downward pump arrangements, 
except that the upward-looking arrangement was easiest to align with 
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the bubble stream. The upward-looking pump was used during most 
data collection because of its easily adjusted components. 

The pump transducer orientation did have a slight effect on 
very slowly rising bubbles. If the bubbles were moving slowly enough, 
the bubble path could be altered very slightly by the pump sound field. 
The altered path was more pronounced for the downward-looking 
pump transducer arrangements. 

The last experimental phase was conducted in the large 
acoustically insulated tank, and resulted in a lower (several dB) back- 
ground noise level across the spectrum. This was significant in that 
the lower resonance frequency (larger radii) sideband pressure levels 
were slightly higher above the noise level than they had been in the 
aquarium. The large tank further demonstrated the need for a reflec- 
tion-free environment to reduce background noise, just as the absorp- 
tive material had reduced reflections and noise in the aquarium. The 
ocean environment and surface ship wakes are quiet in the imaging 
frequency (ultrasonic) region. Use of absorptive material to eliminate 
reflections is also critical to the Dual Frequency Pump Method. 

6. Bubble Counting. Photography, and HP3585 Resolution Band 

Width 

Actual bubble counts were not obtained here using the Dual 
Frequency Pump Method of bubble density measurements for three 
reasons. The first and main reason was that the “max hold” feature on 
the spectrum analj^er, as previously mentioned, was not a method 
sophisticated enough to count bubbles. The second was that the 
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receive transducer was not calibrated and therefore the sideband 
sound pressure amplitudes could not be calculated from the recorded 
sound pressure levels. The third reason was that efforts at photogra- 
phy failed to provide a clear photo enlargement so that a bubble count 
in the bubble stream with a magnifying comparator could be accom- 
plished. Photographs were taken with a Minolta XD-5 35 mm camera 
and a standard 50 mm lens. More sophisticated photographic equip- 
ment, specifically a short focal length lens, is needed to obtain close- 
up photographs. 

Even though actual bubble counts at each sideband frequency 
were not obtained, the plots show the presence of different bubble sizes at 
the discrete coi ± cop locations on the frequency spectrum. The resolu- 
tion bandwidth of the HP3585A spectrum analyzer is critical when deter- 
mining the bubble sizes present using the “max hold” feature. Figure 
set 42 plots the upper dual-frequency sidebands using a 1.0 kHz resolu- 
tion and an extra-wide spectrum window for accuracy. The plot shows 
that small resolution bandwidth provides narrow sound level “spikes” 
that are key to increased precision in bubble size determination. 

Figure 42c shows the presence of discrete bubble sizes rang- 
ing in resonance frequency from 100 kHz to 310 kHz. The pump fre- 
quency cut-off was 310 kHz. Bubble size discrimination is possible, 
but not the actual number of bubbles present. More sophistication in 
photography, sound field amplitude control, transducer types, and 
computer software will provide the means to obtain actual bubble 
counts for each bubble size. 
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B. RECOMMENDATIONS 



The recommendations for further investigation and conduct of 
dual-frequency pump bubble cloud density measurements are listed 
briefly here. The concept is most feasible. With further study and 
development, the design of an acoustical measurement device would 
not be overly difficult. 

1. Only small, narrow bubble streams are practical for measure- 
ment due to bubble screening. The bubble stream dimensions 
must be smaller than that of the sample volume to avoid sound 
field attenuation through bubble screening. Measurement of 
large, dense bubble clouds must be done in a manner which iso- 
lates the sample volume to avoid screening. 

2. The success of dual-frequency bubble cloud measurements rests 
on obtaining the proper signal- to-noise ratio. For this reason, noise 
levels in the sideband frequency range must be kept to a mini- 
mum so that low-amplitude sideband levels can be detected. 

3. The sample volume should be placed in the far field of aU trans- 
ducer beam patterns to prevent the sample volume from being 
in a “null.” 

4. Use statistics as shown in this thesis to determine the duration 
of the sample time. Also, if a large bubble cloud is the target of a 
device, design the device to have several sample volumes and 
statistically average the results. 

5. The problem is complicated enough to warrant the use of a 
computer to run the spectrum anal 3 ^er and plot its data. Inte- 
grate an HP300 series computer into the equipment rack to 
compare and control data. T^is is a most important step which 
is critical to conducting accurate bubble counts. 

6. More sophistication is needed at the pump transducer and 
pump transducer input level. Investigate using broadband noise, 
transducers in series, high-quality flat -frequency response 
transducers, or computer-controlled transducer input power to 
solve the inability to provide constant pump pressure amplitudes 
at the sample volume. 

7. Investigation of higher imaging frequency and higher imaging 
input power is needed to provide better signal-to-noise ratios. 
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Receive transducers with higher center frequencies should be 
used to increase transducer response. 

8. Investigation of the optimum minimum pump pressure ampli- 
tude in the presence of the "widest possible pump frequency 
range (cOp) is needed. The pump frequency (cop) range must be 
as wide as possible to include all bubble sizes at different hydro- 
static pressures. At the same time, the pump sound pressure 
amplitude (pp) must be kept low enough to reduce harmonics 
(which raise the noise level). A balance between pump fre- 
quency cop and amplitude (pp) must be found and at the same 
time provide the bubbles present "with enough sound pressure 
amplitude to produce resonance. 

9. Use of a calibrated receive transducer is critical to determining 
the received pressure amplitudes of the signals from the bub- 
bles. A calibrated receive transducer would allow immediate 
computation of sideband pressure amplitudes and estimates of 
the bubble sizes and numbers present. The computer could also 
perform this function with data from the HP3585A. 

10. Some improved means of visual or photographic verification is 
necessary to confirm the dual-frequency pump method of mea- 
surement. Investigate more sophisticated photography and back- 
ground lighting. Also, investigate other acoustical means (such 
as those developed by Medwin) of bubble density measurements 
to partially validate the Dual Frequency Pump Method. 

In summary, the Dual Frequency Pump Method offers many 
advantages over other acoustic techniques and, when linked to a com- 
puter with the proper software, may be the best technique for mea- 
suring bubble densities in ship wakes. Many hardware improvements 
are needed in the Naval Postgraduate School Dual Frequency effort, but 
the problem lends itself easily to further Acoustics/ Systems Engi- 
neering graduate -level study. Therefore, further investigation and 
study into the Dual Frequency Pump Method is highly recommended. 
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APPENDIX A 



RESONANCE FREQUENCY COMPUTATION PROGRAM 



BUBRAD IS A FORTRAN PROGRAM THAT CALCULATES THE RESONANCE 
FREOUENCY FOR BUBBLES OF A GIVEN RADIUS IN SEAWATER NEAR THE 
SURFACE. 



***************************************************************** it* ********************** ******* 



BUBBLE RADIUS DATA CAN BE FOUND IN FILE: BUBRAD1 DATA A1 
************************************** LI3T OF VARIABLES ♦*♦*♦*♦**♦******♦********♦♦*****♦♦♦♦♦* 



PI = PIE 

GAMMA = RATIO OF SPECIFIC HEATS 
SIGNA = SURFACE TENSION (DYNES/CM) 

ROW = SEAWATER DENSITY (GRAMS/CM3) 

R = MEAN BUBBLE RADIUS 
P = PRESSURE (ATMOSPHERES) 

W1 = NATURAL ANGULAR FREOUENCY 
W2 = NATURAL ANGULAR FREOUENCY (SURFACE TENSION) 
FRE01 = FREOUENCY IN KHZ (WITHOUT SURFACE TENSION) 
FRE02 = FREOUENCY IN KHZ (INCLUDES SURFACE TENSION) 



******************************************************************************************************* 



GAMMA=1 .4 
SIGMA=70.0 
l=0 

ROW=1 .026 
PI=3.141593 
DYNES=1.0133E+6 
WRITE(6,10) 

1 0 FORMAT(/,1X,’INPUT TOTAL HYDROSTATIC PRESSURE, P’) 
READ(5,*) P 
P1=P*DYNES 
WRITE(6,1000) P 
100 READ(1,200,END=9999) R 
200 FORMAT(F7.4) 
l»l+1 

W1 =(1 .0/R)*SORT((3.0*GAMMA*P1 )/ROW) 
FREQ1=(W1/2.0*PI))*0.001 

W2=(1 .0/R)*SORT((3.0*GAMMA*(P1 +(2.0*SIGMA/R))-(2.0*)/R))ROW) 
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FREQ2=(W2/2.0*PI))*0.001 
WRITE(6,2000) R.FREQ1.FREQ2 
IF (I .LT. 10) GO TO 100 

1 000 FORMAT(/,1X,’FOR SEAWATER NEAR THE SURFACE AT PRESSURE P=',1X.F3.1 

*,///, IX, 'RADIUS (CM)', 5X, 'NATURAL FREQ (KHZ)',5X,'FREQ W/SURF TENS 
*ION',/) 

2000 FORM AT(/,4X,F7.4,10X,F9 .4,1 3X,F9.4) 

9999 STOP 
END 
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APPENDIX B 

SPEED AMPLITUDE (Uq) FOR A PULSATING SPHERE 

Recall from Chapter II that the bubble radius R is defined by 
equation (10), which is written 

R = Ro (1 + x). 

Differentiation of equation (10) with respect to time yields 

' dx 
U = R = Ro^. 

where U is the speed of the pulsation. The assumed solution for x is 
written in the form of 



X = Aj cos (cojt + 0). 



Equation (12) expanded on this form in Chapter II to include the 
necessary frequencies of the pump and image sound fields. Now, from 
above, the equation for x is substituted into the equation for U and 
yields 



U = Ro 



dx 

dt 



U = Ro Aj ^ cos (cOjt + 0) 



U = - Ro Aj coj sin (cojt + 0). 
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Pulsation speed U can also be written as a function of the speed 
amplitude U©. Ignoring the negative sign (we are concerned only with 
amplitude), pulsation speed is written 

U = Uo sin (cojt + 0). 

Therefore, the speed amplitude is 



Uo = Ro CDj. 
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APPENDIX C 



UPPER SIDEBAND PRESURE AMPLITUDE 
COMPUTATION PROGRAM 



C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



SIDEBND IS A FORTRAN PROGRAM FOR COMPUTING THE UPPER DUAL FREQUENCY 
SIDEBAND PRESSURE AMPLITUDES (DYNES/CM2) FOR BUBBLES RADIATED BY 
THE DUAL FREQUENCIES USING EQUATIONS (28A) AND (28B) AT 1.0 ATMOSPHERES. 



BUBBLE RADIUS DATA CAN BE FOUND IN FILE: BUBRAD1 DATA A1 

list of variables 






RO = MEAN BUBBLE RADIUS 

U = VISCOSITY 

ROW = SEAWATER DENSITY 

W1 - NATURAL ANGULAR BUBBLE FREQUENCY 

W1 = IMAGING FREQUENCY 

WP = PUMP FREQUENCY 

PI = IMAGE PRESSURE (DYNES/CM2) 

PP = PUMP PRESSURE (DYNES/CM2) 

R1 = DISTANCE TO THE BUBBLE (CM) 

DELTA = DAMPING COEFFICIENT 

X2 = DIMENSIONLESS FREQUENCY AND DAMPING TERM 
PPL1 = SUM SIDEBAND PRESSURE AT RESONANCE (DYNES/CM2) 
PPL2 = SUM SIDEBAND PRESSURE NOISE (DYNES/CM2) 






U=0.01 

R1=6.7 

1=0 

GAMMA=1.4 

ROW= 1.026 

DYNES=1.0133E+6 

PP=830.8 

Pl=659.9 

WRITE(6,1000) 

100 READ(1,200,END=9999) RO 
200 FORMAT(F7.4) 

1 = 1+1 

W1=(1.0/R0)*SQRT((3.0*GAMMA*DYNES)/ROW) 
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F1=(W1/(2.0*3.141593))*0.001 
DELTA={4.0*U)/(ROW* W1 *(R0**2)) 

WP=W1+W1/50.0 

0HM2=WP/W1 

Z1=(ROW*W1*RO)**2 

X2=1 ,0/SQRT((1 .0-OHM2**2)**2+(DELTA“2)*(OHM2**2)) 

PPL1=(R0W*PI*PP*R0)/(Z1*DELTA*R1) 

PPL2=(ROW*PI*PP*RO*X2)/(Z1*R1) 

WRITE(6,2000) R0,F1,PPL1,PPL2 
IF (I .LT.10) GO TO 100 

1000 FORMAT(/,1X.'RADIUS (CM)',5X;SUM PRESSURE • PEAK’.SX.’SUM 
PRESSURE NOISEV) 

2000 FORMAT(/,2X,F7.4.12X.F9.4.10X.F9.4.14X.F11.8) 

9999 STOP 
END 
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APPENDIX D 



PUMP TRANSDUCER PRESSURE AMPLITUDE CALCULATION 

In order to calculate the sound pressure amplitudes supplied by 
the pump transducer and measured at the sample volume by the LC-10 
probe transducer, it is necessary to formulate two conversion formulas. 
The first conversion is for the plots from the HP3585A, shown in 
Chapter IV. The HP3585A spectrum analyzer plotted transducer 
response in units of dBm (dB re 1 mW), therefore, the first conversion 
changes the plotted output signal of the LC-10 from dBm to volts. 
. Because dBm is referenced to milliwatts, 

p 

dBm = 10 logic 

vdiere Pref =1.0 mW. Substituting P = V^/R yields 
dBm = 10 logio (^] = 10 logic 

Rearranging the terms gives the equation which converts dBm to volts: 

V 

,V(Pref)(R)J ■ 

where 



dBm = 20 logic 
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Two examples of these useful calculations would be to calculate 
the sound pressure from the pump transducer at the location of 
sample volume approximately 10 cm from the transducer face. From 
Figure D-2"at 100 KHz, the LC-10 sensed approximately -51 dBm at 
10 cm distance away , from the pump transducer using the HP3585A 
with the reference resistance (R) of 50 Q. The pump sound field 
sweep time used to plot Figure D-3 was 0.01 seconds. Now, 

- 43.0 dB„ = 20 logio [; ?a.Tx T5:3w)'('50n). 

)delds V = 1.583 x lO”^ volts. Therefore, pump sound pressure pp 
measured at the sample volume is 

Pp = VlAf = (1.583 X 10-3 volts) |^1.4125 x 105 



Pp = 223.6 dynes/cm^. 

If the pump sound field is swept every 0.10 seconds, the sound 
pressure level increases, in this case by 11,4 dBm. From Figure D-4, 
at 100 KHz the equation 

- 31.6 dBm = 20 logio [^(i.oxlO-3) (50 £l). 

yields v = 5.88 x 10“3 volts. Therefore, 



Pp - 830.8 dynes/cm^. 
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Figure D-2. Pump Transducer Response at 10.0 Centimeters— Pump Sweep 0.01 Seconds 
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Figure D-3. Pump Transducer Response at 10.0 Centimeters — Pump Sweep 0.10 Seconds 
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Figure D-4. Pump Transducer Response at 7.0 Centimeters— Pump Sweep 0.10 Seconds 



The same set of equations is used to calculate the pump pressure 
amplitude at 7.0 cm from the pump transducer face. From Figure D-4, 
with a pump frequency sweep time of 0.10 seconds, the sound 
pressure level at 100 KHz is - 39.2 dBm. Therefore, pp = 346.3 

d3oies/cm2. 
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APPENDIX E 



IMAGE TRANSDUCER PRESSURE AMPLITUDE CALCULATION 

Appendix D calculated the pump transducer’s pressure amplitude 
(Pp) at the sample volume using the sound pressure level recorded on 
the HP3585A and knowing the calibrated frequency response of the 
LC-10 hydrophone. The image transducer’s pressure amplitude (pi) 
at the sample volume is calculated in the same fashion, only here the 
sound pressure level is measured with an ultrasonic hydrophone, the 
MA-1, designed by Bob Bruce and Bob Middleton of the Naval 
Postgraduate School, Monterey, California. The ultrasonic hydrophone 
calibrated frequency response is shown in Figure E-1. The sensitivity 
levels of 2.25 MHz and 2.50 MHz are -97.0 dB and -97.5 dB, 
respectively. Using the equation from Appendix D, 

ML = 20 logio ^ . 

5Wref 

and that the reference level Mxti is 1.0 volts/Pa, the conversion factor 
M is calculated as 

M (2.25 MHz) = 7.0795 x lO^ ^ = 7.0795 x 105 



and 



M (2.50 MHz) = 7.4989 x lO^ ^ = 7.4989 x lO^ ^T^ . ^s/cm^ 
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Figure E>-1. MA-1 Hydrophone Pressure Sensitivity 



When the imaging transducer was driven at 6.0 volts, the ultrasonic 
hydrophone sensed a level of -47.6 dBm and -48.1 dBm (using the 
HP3585A Spectrum Analyzer) at 2.25 MHz and 2.50 MHz, 
respectively, in the imaging far field sample volume (see Figure E-2). 
Now, at 2.25 MHz, 



- 47.6 dBm = 20 logic 



V 

,V(1.0x 10"^w)(50Q) 



)delds V = 9.3215 x lO'"*^ volts. Therefore, imaging sound pressure pi 
at the far field sample volume is 

Pi = VfW = (9.3215 X 10-4 ) f7.6795 x 10^ ^y^^s/cm^ 



Pi = 659.9 dynes/cm^. 

Likewise, the imaging sound pressure pi at 2.50 MHz is 
Pi = 659.9 d3mes/cm^. 

Slightly different transducer frequency responses and hydrophone 
sensitivities produced the same sound pressure. 
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2.150 2.200 2.250 2.300 2.350 2.400 2.450 2.500 2.550 2.600 2.650 

Frequency (MHz) 

Figure E-2. Imaging Transducer Response at 6.7 Centimeters 



APPENDIX F 



CALCULATION OF INSONIFIED SAMPLE VOLUME 




Figure F-1 

Sample Volume Definition 

The sample volume formed by the intersection of the main lobes 
of the transducer beams (3 dB down) can be considered an approxi- 
mate cylinder of radius a or b and length 2a or 2b. Table II shows for 
the 2.25 MHz transducers 

e = 3.10® r = 4.74 cm 

where r is the distance to the center of the sample volume. 

Therefore, a « b « 4.74 tan(3.10°) = .257 cm. Now, volume = area 
times length = (7ta2)(2b). If a * b, then Volume = = 2 k (.257 cm)3 

= .107 cm3. Sample volume is approximately .107 cm3. 
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APPENDIX G 



STATISTICAL SAMPLE TIME 



Calculation of the sample time or duration of the “look” begins 
with knowing how many samples are needed. To obtain a 95-percent 
confidence factor with a maximum error (E) of 5 percent, use 
E = .05 

Za/2 ~ 1.96 

from the normal distribution tables. Equation (41) is now used and 
5 Tields 



The number of samples, n = 384.2, is rounded up to n = 385. 

From Appendix 5, the length of the sample volume is 2b = 
.514 cm. A length must be sampled 385 times for a total sample 
length of 

sample length = (n) (2b) = 197.9 cm. 

Bubbles between 50 and 100 micrometers (pm) rise at a rate of .75 
and 2.1 centimeters /second, respectively [Ref. 23]. If the bubbles are 
rising, for example, at a rate of 1.5 cm/s through the sample volume, 
then 

Total Sample time = 131.9 seconds. 
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Smaller bubbles rise even more slowly, so an even greater sample time 
would then be required to get the. statistically correct number of 
samples. 
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